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PREFACE 


Stanford  Research  Institute  has  been  conducting  a  study  of  position 
location  and  navigation  for  the  U.S.  Marine  Corps  under  Contract  Nonr  2332(00) 
An  important  part  of  this  study  is  the  assessment  of  present  and  future 
systems.  In  developing  criteria  for  such  assessment,  considerable  diffi¬ 
culty  was  experienced  in  locating  suitable  references  for  the  necessary 
mathematical  analyses.  Further,  descriptions  of  systems  accuracy  found 
in  the  literature  were  found  to  be  often  confusing,  sometimes  ambiguous, 
and  occasionally  in  error.  Thus  an  extensive  effort  became  necessary  to 
collect  suitable  information,  to  develop  additional  methods  of  analysis, 
and  to  select  a  uniform  method  of  specifying  system  accuracy. 

This  Heseurch  Memorandum  has  been  prepared  to  record  these  mu thema t ic a  1 
techniques  us  they  have  been  developed  and  used  for  the  assessment  of  var¬ 
ious  position  location  and  navigation  systems.  A  concomitant  result  hus 
been  the  spec i fi c a t i on  a f  useful  methods  of  describing  the  accuracy  of 
position  measurements.  This  memorandum  follows  a  format  of  presenting 
the  results  with  illustrative  examples.  These  are  followed  by  extended 
discussion  and  derivation  of  formulas  in  the  several  appendixes.  The 
final  report  of  the  position  location  and  navigation  study  is  classified. 

Since  the  mathematical  techniques  developed  for  this  study  can  be  applied 
to  other  types  of  problems,  they  have  been  presented  in  a  separate  unclas¬ 
sified  volume  to  permit  wider  dissemination  of  these  techniques  to  those 
who  may  find  them  useful. 

Additional  volumes  dealing  with  this  subject  may  be  issued  under  this 
and  other  studies  as  the  need  urises. 

Typographical  errors  existing  in  the  first  printing  of  this  memorandum 
have  been  corrected  in  this  printing.  No  other  textual  changes  have  been 
made  in  this  second  printing. 
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LOP  Line  of  position 

At  u  point  determined  by  the  intersection  of  two  lines  of  position 

c  Standard  deviation  associated  with  LOP  #1,  measured 
perpendicular  to  the  LOP 

cr  Standard  deviation  associated  with  LOP  #2,  measured 
perpendicular  to  the  LOP 

d  Angle  between  the  two  lines  of  position,  also  (and 
equi vulently) angle  between  LOP  1  and  LOP  2  as  shown 
be  low. 


cr*  Fictitious  standard  deviation 

d*  Fictitious  intersection  angle 

The  combination  of  cr*,  cr* ,  and  a*  represents  an 
equivalent  description  of  a  probability  distribu¬ 
tion  actually  described  by  cr^,  o-  and  a. 

After  transformation  of  a  probability  distribution  described  by  Cj,  cr2 , 
and  a  to  an  equivalent  distrioution  defined  in  terms  of  the  axes  of  an 
e  1  i  ipse 

< y standard  deviation  along  major  axis 
o  standard  deviation  along  minor  axis 
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When  distributions  ubout  several  points  are  combined  to  obtain  un  overull 
distribution,  double  subscripts  ure  used  to  designate  the  separate  ellipses 
which  describe  the  distributions  ubout  each  point. 

Standard  deviation  LOP  // 1 — ellipse  ff  1 
Standard  deviation  LOP  1} 2 — ellipse  II  1 
crl  2  Standard  deviation  LOP  — ellipse  //  2,  etc, 

After  transformation  to  standard  deviations  along  the  axes  of  the  indivi¬ 
dual  ellipses. 


Standard 

deviation 

along  major 

axis— el  1  ipse 

in 

V 

Standard 

deviation 

along  minor 

axis — e 1 1  ipse 

tn 

at  2 

Standard 

deviation 

along  major 

ax i s — e 1 1  ipse 

in, 

*1 

Angle  between  x-axis  of  ellipse 
coordinate  axes  for  combination 

H 1  and  ar hi t 
of  ell ipses 

r  ar  y 

1 ’  aii  Standard  deviations  along  the  arbitrarily  selected 
a*  2  ,  o  ^  2  axes,  designated  ui  and  2,  for  multiple  ellipse 
etc.  combination, 

Pj  ,  p  j,  Cross  product  function  involved  with  trails  formation 
etc.  to  w  and  z  axes. 


In  general  discussion  of  severul  ellipses  the  subscript  i  is  used  to 
designate  the  general  function.  /  »  1,2,3 . 

In  the  final  combined  ellipse,  the  subscript  f  is  used. 
o  As  above  for  the  final  ellipse 


cr 


«  / 


Standard  deviation  along  major  axis 
Standard  deviation  along  minor  axis 
Angle  between  x-axis  and  arbitrary  coordinate  axes. 


1.  BASIC  STATEMENT  OF  POSITION  LOCATION  PROBLEM 


The  busic  problem  in  position  locution  is  the  determination  n  i'  the 
coordinates  of  u  remote  point  with  respect  to  u  known  or  arbitrury  refer- 
enee.  The  remote  point  may  be  a  landing  /.one  for  troops,  or  a  target  upon 
which  it  is  desired  to  deliver  ordnance.  Many  other  examples  will  immedi¬ 
ately  come  to  mind,  The  problem  giving  rise  to  the  necessity  for  complex 
mathematical  unalysis  is  the  fact  that  no  measurement  can  be  made  without 
error.  Thus,  the  results  of  a  position  determination  in  fact  must  be 
described  in  terms  of  the  probability  of  being  within  a  given  distance  of 
the  desired  point.  Actually,  this  last  statement  is  in  too  simple  terms: 
because  more  than  one  error  is  usually  involved  in  the  sum  total  of  meos- 
inetits,  it  becomes  necessary  to  consider  the  shape  of  the  probability  dis¬ 
tribution  about  the  desired  point.  In  general,  these  probability 
distributions  are  ellipses  rather  titan  circles. 

A  relatively  simple  problem  in  position  location  is  given  in  Fig.  1. 
Here  is  assumed  a  reference  baseline  established  by  the  measurement  system, 
The  location  of  an  artillery  battery  is  then  measured  by  the  system,  giving 
rise  to  an  ellipse  within  which  a  given  probability  muy  be  stated  that  the 
artillery  buttery  is  actually  locuted.  Then,  with  respect  to  the  same  base 
line,  u  forward  observer  is  located,  giving  rise  to  a  second  ellipse  within 
which  he  may  be  located  to  a  stated  percentuge  probability.  From  his  loca¬ 
tion  he  makes  measurements  on  a  target  which  then  may  be  located  within  a 
still  different  ellipse.  From  this  information  the  dotted  line  giving 
firing  orders  in  range  und  azimuth  is  calculated.  The  weapon  effects 
ellipse  is  shown  dotted  and  superimposed  upon  the  target  location  ellipse. 
Then,  the  problem  of  immediate  interest  is  to  calculate  the  probab i  1 i ty o f 
damage  to  the  target.  Techniques  have  been  established  to  perform  such 
calculation  when  each  of  the  error  figures  about  the  various  points  in  the 
problem  is  a  circle.  However,  the  use  of  such  circles  can  be  quite  mis¬ 
leading  when  the  actual  figures  are  ellipses.  In  particular,  the  weupons 
effect  pattern  is  commonly  a  very  elongated  ellipse,  differing  greatly 
from  a  circle.  It  is  also  characteristic  of  many  electronic  measuring 
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FIG.  1  BASIC  LOCATION  PROBLEM 


techniques  that  the  results,  to  be  meaningful,  rr.usL  bo  expressed  in  terms 
of  ellipses.  Thus  it  becomes  necessary  to  develop  u  procedure  of  analysis 
which  permits  the  consideration  of  ellipses. 

A  circle  is  reudily  defined  in  terms  of  its  center  and  a  single 
distance,  the  radius.  An  ellipse  requires  udditionul  information  —  the 
center,  and  two  distances  to  correspond  to  the  radius  of  a  circle — the 
semi-major  axis  and  the  semi-minor  axis.  Further,  when  we  are  concerned 
with  more  than  one  ellipse,  we  must  also  be  concerned  with  their  relative 
orientations.  Thus,  any  analytical  procedure  concerned  with  the  end  result 
of  the  consideration  of  a  number  of  error  ellipses  considers  the  elements 
of  each  ellipse  and  the  angles  of  these  elements  with  respect  to  a  common 
reference. 

Because  the  detailed  consideration  of  analytical  techniques  concerned 
with  any  single  ellipse  is  itself  quite  complex,  the  discussion  will  detail 
first  the  considerations  of  a  single  ellipse  of  error.  Following  this 
presentation,  the  method  of  combining  several  ellipses  to  obtain  the  end 
result  will  be  described. 
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2.  ANALYSIS  OF  A  SINGLE  ERROR  ELLIPSE 


Most  oi'  the  position  locution  systems  considered  in  this  program 
determine  the  locution  of  u  point  ut  the  intersection  of  two  lines  of 
position.  However,  both  lines  of  position  may  be  in  error.  Figure  2 
shows  such  an  intersection  of  two  lines  of  position.  The  lines  of  po¬ 
sition  in  this  illustration  are  range  measurements  from  two  points  at 
the  extremities  of  a  baseline  of  known  length.  The  measurements  of 
ranges  are  reported  as  some  numerical  value.  However,  because  of  inac¬ 
curacies  in  measurement,  the  actual  range  may  not  be  the  indicated  value, 
but  may  lie  somewhere  between  the  limits  shown  as  additional  arcs  either 
side  of  the  measured  line,  Thus,  one  becomes  interested  in  the  proba¬ 
bility  that  the  uciuul  point  lies  within  some  close  distance  of  the 
indicated  point. 

The  intersection  of  the  two  lines  of  position  together  with  the 
standard  deviations  associated  with  each  are  shown  to  expanded  scale  in 
Fig.  3.  (Standard  deviation  is  a  measure  of  error  and  is  defined  in 
Appendix  A.)  Standard  deviation  as  a  measure  of  error  is  commonly  desig¬ 
nated  by  the  Greek  letter  sigma  ( cr )  and  the  Greek  letter  alpha  (ft)  will 
be  used  throughout  the  analysis  to  designate  the  angle  of  intersection  of 
two  lines  of  position.  It  can  be  shown  that  the  contours  of  equal  prob¬ 
ability  density  about  such  an  intersection  are  ellipses  centered  about 
the  intersection  of  the  two  lines  of  position.  Thus,  the  ellipse  shown 
in  Fig.  3  might  be  the  75%  probability  ellipse,  meaning  that  there  is  a 
75%  probability — three  chances  in  four— that  the  actual  position  of  the 
point  whose  location  is  desired  lies  within  the  ellipse  drawn. 

The  detuiled  statistical  analysis  of  the  diagram  shown  in  this  figure 
is  quite  complex.  Salient  features  of  the  analysis  will  be  stated  in  the 
main  portion  of  this  memorandum.  Higorous  mathematical  analyses  will  be 
found  in  the  appendixes.  This  main  body  of  the  memorandum  will  indicate 
the  assumptions  made  in  the  analysis  of  the  accuracy  of  position  location 
systems  and  will  present  formulas,  graphs,  and  nomograms  to  obtain  numer¬ 
ical  results.  Typical  examples  will  be  illustrated. 
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FIG.  2  POSITION  LOCATION  AT  INTERSECTION  OF  TWO  LINES  OF  POSITION 
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FIG.  3  EXPANDED  VIEW  OF  INTERSECTION  OF  TWO  LINES  OF  POSITION 


In  t  lie  unulysis  In  follow  of  the  uerurucy  of  position  locution  systems 
tlic  following  a  s  s  mu  p  t.  i  on  s  have  boon  mu  do  : 


(1)  All  bius  errors  hove  been  removed,  leuving  only  the  rundom 
errors  to  be  anulyzed.  In  muthemuticul  terms,  the  mean  or 
uverugc  error  is  ussumed  to  be  zero. 

(2)  These  rundom  errors  are  assumed  to  be  normally  distributed. 

The  mat  hema i  e  a  1  implications  of  this  assumption  are  dis¬ 
cussed  in  following  paragraphs.  This  assumption  is  required 
to  permit  a  mu  t  lienia  t  i  cu  1  .statement  of  overull  accuracy 
relationships  to  be  developed. 

(3)  The  errors  associated  with  the  two  lines  of  position  are 
assumed  to  be  independent.  This  assumption  implies  that 

a  change  in  the  error  of  one  lint:  ol'  position  has  no  effect 
upon  the  other.  This  assumption  permits  the  analysis  to 
consider  unequal  errors  associated  with  the  two  lines  of 
position,  thus  assuring  a  realistic  .mathematical  model. 

(4)  The  lines  of  position  are  assumed  to  be  straight  lines  over 

the  small  area  of  interest  in  the  lie  i  gli  bur  hood  of  the  point,  the 
position  of  which  is  desired.  (See  Fig.  3.)  This  assump¬ 
tion  is  valid  so  long  us  the  stundard  deviation  is  small 
with  respect  to  the  uetuul  radius  of  curvature  of  the  line 
of  position.  In  the  unulysis  of  systems  this  is  usually 
found  to  be  the  cuse.  Not  to  muke  this  assumption  would 
unreasonably  complicate  the  mathematical  unulysis. 

(5)  The  unulysis  of  errors  of  position  is  confined  to  two  di¬ 
mensions.  The  third  dimension,  ultitudt,  may  be  considered 
separately,  if  desired,  if  the  system  being  analyzed  is 
capable  of  altitude  measurements. 

As  shown  in  Fig.  3,  the  general  case  of  the  intersection  of  two  lines 
of  position  at  any  angle  and  with  different  values  of  error  associated 
with  each  line  of  position  results  in  an  elliptical  error  figure.  Simpli¬ 
fied  to  geometrical  terms,  the  ellipse  looks  like  that  of  Fig.  4. 

From  this  illustration,  one  may  readily  surmise  that  the  exact  shupe 
of  the  error  figure  varies  with  the  magnitudes  of  the  two  input  errors, 
o j  und  cij,  as  well  as  with  the  intersection  angle,  (i.  The  unglc  Ot  is  also 
the  angle  between  the  two  values  of  sigma  because  the  standard  deviations 
are  mutually  perpendicular  to  their  corresponding  lines  of  position.  How 
these  variations  may  be  calculated  was  the  objective  of  considerable  ana¬ 
lytical  effort  by  the  project  team  after  a  literature  survey  indicated  the 
inadequacy  of  available  techniques. 


8 


Results  of  tliis  unulyticul  effort  have  given  two  separate  but  inter¬ 
related  methods  of  obtaining  the  desired  results.  Each  will  be  described 
us  simply  us  possible  und  the  necessary  formulas  given  with  examples, 
lletuiled  deviations  of  the  formulas  obtained  may  be  found  in  the  appendixes. 

Using  either  of  the  two  methods  of  anulysis  developed,  the  end  result 
is  the  determination  of  the  probability  that  the  point  is  Located  within  a 
circle  of  stated  radius.  The  basis  of  this  concept  may  best  be  seen  by 
considering  for  u  moment  the  special  case  when  the  two  errors  are  equal 
and  the  utigle  of  intersection  of  the  lines  of  position  is  a  right  angle. 

In  this  case,  and  in  this  cuse  alone,  the  error  figure  becomes  a  circle 
and  is  described  by  the  circular  normul  distribution.  A  plot  of  this 
special  function  is  given  in  Rig.  5.  The  plot  is  to  be  interpreted  as 
follows:  the  horizontal  axis  is  measured  in  terms  of  fl/cr,  fl  being  the 

radius  of  a  circle  within  which  it  is  desired  to  be  located,  and  a  being 
the  error  measure,  The  error  measure  is  given  simply  as  cr ,  for  in  this 
circular  case  cr  t  =  o j.  To  illustrate,  a  measurement  system  gives  a  cir¬ 
cular  error  figure  and  has  u  vulue  of  cr  «  10  meters;  the  probability  of 
actually  being  located  within  a  circle  of  10  meters  radius  when  fl/cr  »  1.0 
may  be  read  from  the  vertical  uxis  to  be  39.3%.  To  obtain  the  CEP,  the 
radius  of  a  circle  within  which  a  50%  probability  results,  the  corre¬ 
sponding  value  of  R/a  is  seen  to  be  1.18  from  the  graph.  Thus,  for  this 
example,  the  CEP  would  be  11,8  m.  The  concept  of  the  function  fl/cr  will 
be  found  useful  in  more  complex  cases, 
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FIG.  5  CIRCULAR  NORMAL  DISTRIBUTION 
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Slurting  with  the  inputs  of  Fig.  4,  it  is  assumed  that  it  is  possible 
to  find  fictitious  values  of  sigma  so  that  the  two  differing  values  origi¬ 
nally  given  may  be  replaced  by  two  new  sigmas  of  identical  value,  indi¬ 
cated  as  <J* .  At  the  same  time  u  new  and  fictitious  angle  of  intersection 
0( *  is  also  required.  Figure  6  indicates  these  new  values. 


FIG.  6  TRANSFORMED  PARAMETERS  OF  ERROR  ELLIPSE 


Also  required  for  use  with  this  first  method  is  a  whole  set  of  prob¬ 
ability  curves,  similar  to  that  of  Fig.  5,  but  with  a  separate  curve  for 

each  value  of  intersection  angle,  Such  curves  have  been  calculated  with 

the  aid  of  a  digital  computer  and  are  shown  on  Fig.  7.  These  curves  can 
be  used  only  when  the  two  error  measures  urc  equal,  hence  the  need  for 
making  the  transformation  of  the  previous  paragraph. 

The  vulues  of  the  functions  cr*  and  a*  needed  to  utilize  the  curves 
of  Fig.  7  may  either  be  computed  from  the  formulas  given  below  or  may  be 
determined  from  Figs.  8  and  9, 

s  i  u  2/Vcr*  +  oj 

O'*  =  • 

a*  =  a  resin  (sin  2/3  sin  «) 


wli*.'  ri‘ 


(i  *  p rctan  (o-j/o-j) 

t  Ihi  ft 

2crlo-2 

a  i  n  2  fi  =  - — 

■d  *  <,> 

Tlit*  do  r  i  vu  t  i  on  of  these  formulas  may  be  found  in  Appendix  D.  It  is  ob¬ 
vious  I. hut  considerable  computation  is  required  to  use  the  formulas 
directly,  hence  the  curve  of  Fig.  8  and  the  nomogram  of  Fig.  9  will  be 
found  to  facilitate  the  use  of  the  formulas. 

First  one  must  calculate  the  ratio  ctj/ctj.  cr  l  is  always  taken  as 
the  larger  of  the  two  in  this  fraction,  such  that  the  value  is  always 
less  than  1.0.  With  this  ratio,  enter  the  curve  of  Fig.  8  and  obtain  the 
cr*  factor.  Multiply  &  by  this  factor  to  obtain  the  fictitious  function 
cr* ,  The  nomogram  of  '  9  is  used  with  the  sane  ratio  to  obtain  the 

fictitious  ungle  Ot*. 

A  numerical  example  will  illustrate  the  method  of  calculation. 

Assume  a  position  location  system  has  provided  the  following  data: 

a  »  50° 

v j  «■  20  meters 

Oj  *  15  meters 

Whut  is  the  probability  of  the  location  of  the  point  within  a  circle 
of  30  meters  radius? 

Culculute  the  ratio  cr  2/a  j  *  15/20  ”  0.75. 

Enter  the  curve  of  Fig.  8  with  this  value  and  obtain  the  cr* 
factor  B  0.845.  Multiply  this  value  by  ctj  *  20  to  obtain  cr*  -  16.9  meters, 
Calculate  the  ratio  R/c*  -  30/16.9  ■  1.78. 

Enter  the  nomogram  of  Fig.  9  with  the  same  ratio  cr 2 /cr j  =  0.75  and  with 
the  given  angle  (X  *  50°  to  obtain  the  fictitious  angle  ct*  *>  47°. 

The  values  fl/cr*  «  1.78  and  a*  »  47°  may  be  entered  on  the  curves  of 
Fig.  7  to  obtain  P  •  0.62  or  62%,  interpolating  visually  between  the  40° 
and  50°  curves  for  a*  -  47°. 
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Checking  the  values  obtained  from  the  nomograms  by  substituting  the 
figures  into  the  exact  formulas  we  obtain: 


sin  2/3 


2  x  15  x  20 
225  +  400 


600 

625 


0.96 


iin  50°  a  0.766 


-  u resin  (0.96  x  0.766) 

*  arcsin  0.735 

«  47.3°  (compared  with  47°  from 

the  graphical  solution) 

0.96/623 

/? 

-  0.707  x  0.96  x  25 

-  16.96 


f?/<7*  »  30/16.96  ~  1.77  (compared  with  1.8 

from  the  graphical 
solution) 


The  values  obtained  from  the  figures  are  seen  to  be  sufficiently 
accurate  for  comparative  evaluation. 

Because  the  curves  of  Fig.  7  tend  t.o  crowd  quite  closely  together 
for  some  values  when  the  angles  approach  90  degrees,  Table  I  presents 
the  same  data  in  numerical  form. 


2.2  METHOD  2 

A  second  method  of  working  with  the  error  ellipses  starts  with  a 
different  transformation.  In  this  method  new  values  of  sigma  are  found 
along  the  major  and  minor  axes  of  the  ellipse  according  to  the  formulas 
given  below.  The  geometrical  considerations  are  given  in  Fig,  10. 
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P(R/<r,a.)  CIRCULAR  ERROR  PROBABILITIES 
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FIG.  10  TRANSFORfMT ION  TO  STANDARD  DEVIATIONS  ALONG  ELLIPSE  AXES 

+  o'*  +  +  cr^ )  -  4  s i ii 2  a  rr  Jo- * 

+  oj  -  i /(cr*  -r~a22)  “  4  sin2  acr^o^ 

The  complex  derivation  of  these  formulas  is  given  in  its  entirety  in 
Appendix  B.  Note  that  these  formulas  are  given  terms  of  variances — 
squares  of  the  standard  deviations. 

After  the  values  of  sigma  along  the  orthogonal  axes  of  the  ellipse 
have  been  obtained,  the  re.ults  of  computations  obtained  by  Harter  (Ref.  4) 
may  be  utilized  to  obtain  desired  circles  of  probability.  To  utilize 
Harter's  data  it  is  first  necessary  to  compute  the  ratio  c  «  cr^/a*  where 
cr x  is  the  larger  of  the  two  new  standard  deviations  just  computed.  Harter 
presents  tables  (Tables  II  and  III)  which  then  relate  ellipses  of  varying 
values  of  ellipticity  to  the  radii  of  circles  of  equivalent  probability. 

The  use  of  this  second  method  will  be  shown  by  using  the  same  example 
given  as  an  illustration  for  the  first  method. 
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.1742045 

.1451808 

.1237982 

.1076237 

.0050495 

.0850320 

0768837 

0.3 

1  .3820240 

.3756884 

.3481700 

,3003001 

.2532953 

.2152880 

.1857448 

.  1020820 

.  1443941 

.1206280 

.1175031 

o.o 

.4614038 

.4467708 

I  .4255005 

.3846374 

.3357384 

.2014082 

.2548177 

.2251114 

.  2000797 

.1811783 

1047298 

l) .  7 

j  .5100727 

.5115048 

.4900683 

.4033268 

.4170802 

.3009305 

.3280302 

. 2025064 

. 2029373 

.2381683 

2172055 

0.8 

1  .6702602 

6726967 

.5604457 

,6349387 

.4941882 

.4474207 

.4026628 

.3027122 

.3283463 

.2080700 

.2738510 

0.0 

!  .6318707 

.6288721 

.6191354 

,6093140 

.5851504 

.5213998 

.4759375 

. 4333628 

.3053270 

.3620135 

3330232 

1.0 

.0830808 

.6802326 

,6723580 

.6508242 

,0291249 

.5900953 

.5461319 

.5025790 

.4021421 

.4257553 

.3934093 

1.1 

.7280070 

.7200697 

.7202082 

.7070081 

.0860367 

.0524489 

.6116310 

.5687407 

.6272402 

.4887873 

.4539250 

1.2 

.  7098007 

.7682216 

.7030306 

.7632175 

,7369568 

.7079073 

.0714200 

.6300)08 

.  0803494 

.5408730 

.5132477 

1.3 

.8003990 

.8060648 

,8008554 

,7029968 

,7703650 

.7587266 

.7240073 

.0873122 

.6474304 

.0070822 

.5704426 

1.4 

.8384807 

.8374049 

.8340018 

,8277048 

.3160851 

,7980288 

.7720880 

,738308!) 

.7007000 

.6023035 

.0246889 

1.6 

.8003860 

.8055127 

.8027728 

,8577882 

.8493071 

,8350816 

.8129287 

.7833002 

.7480600 

.7122640 

.6763476 

1.0 

.8004014 

.8897008 

.8875000 

.8834914 

.8708844 

.8667559 

. 8478393 

. 8220240 

.7917104 

.7574708 

.7219027 

1.7 

.0108001 

.0103102 

.908661!) 

,0053700 

.0001740 

,8916630 

.8773110 

.8602471 

.8201137 

.70(7882 

.7402639 

.8021013 

1  .8 

.0281304 

.9270904 

.0263126 

.9237989 

,0197275 

,0130680 

,0010110 

. 8840024 

.8613238 

.8332176 

1  9 

.0426000 

.0422182 

1  .9411290 

.9301680 

.9369865 

.9303018 

.9222277 

1  .9083600 

.8886731 

.8030140 

.8365256 

2.0 

.0544007 

.95422/2 

.0633775 

,961  8416 

.9493015 

.9464540  1 

1 

.0388418 

. 0278700 

.0115702 

.8C014U6 

.8040047 

2.1 

.9042712 

.9040608 

.0034011 

.9622127 

,9603170 

.9573206 

.9522999 

.0437608 

.9305013 

.0122714 

.8807495 

2.2 

1  .9721031 

.9720304 

.9716237 

.9708109 

.9691597 

.9608846 

.9031017 

.0606522 

.0450380 

.9300821 

.9110784 

2.3 

i  .0785618 

.0784275 

.0780408 

.9773460 

.9762419 

.9745230 

.0716034 

,060/300 

.9583739 

.9458085 

.0280946 

2.4 

.9838049 

!  .9835108 

.0832180 

.9820918 

.9818594 

.9805703 

, 97B4001 

.0747496 

.9082698 

.0580804 

.0438652 

2.6 

.9876807 

.9875100 

.9872900 

.9808963 

.9802730 

.9853112 

.9837509 

.0810036 

.9700522 

.9679136 

.9560631 

2.0 

.9000770 

.9000249 

.9904012 

.0901074 

.9897046 

.0889934 

.0878527 

,0858331 

.9821023 

.0756989 

.9660526 

2.7 

1  .9030061 

.9030271 

.9929002 

.0920804 

,0023183 

.0918260 

.0009044 

.0805208 

.0807530 

.9817837 

.9738786 

2.8 

.9948807 

.9948012 

.9047727 

.0040141 

.0943049 

.9939842 

.9933821 

.9923240 

.1)902888 

.9864870 

.9801589 

2.0 

.9002084 

.9002477 

.9901834 

.9080684 

, 9963878 

.9960126 

.9951798 

,9944218 

.9029452 

,9900803 

.9850702 

3.0 

.9973002 

.9072863 

.9972301 

.9071564 

.9970200 

. 3968204 

.9935206 

.0959854 

.0949774 

.9927926 

.0888910 

3.1 

.9980018 

.9080642  | 

.9980212 

.9979022 

.0978099  | 

.9977290 

.9975109 

.9971348 

.9903861 

.0948168 

.9018113 

3.2 

.9086267  | 

.9080182 

.9985949 

.9985633  I 

.9984880 

.9983892 

.0982356 

.0079733 

.0974478 

.9963105 

.9940240 

3.3 

.0090332 

.9090270 

.9990110 

.9089824 

.9939368 

1  .0988077 

.9987607 

j  .9985792 

.9982147 

.0974004 

.9066822 

3.4 

.0993261 

.9093225 

,  .9993112 

.9992909 

9092593 

1  .9992115 

!  .9091376 

.9990120 

.9087620 

.0981808 

.9000113 

3.5  | 

.0996347 

.0096323 

.9995245 

.0995106 

.0994888 

.9994559 

.0094053 

.9093204 

.9991502 

.9987480 

.9978125 

3  0 

.  1)1)903 1 8 

.9090801  | 

.9996718 

.9906053 

.B9905Q5 

.9990281 

.0995938 

.9995304 

.9994218 

.9991442 

.9984002 

3.7 

.0007844 

.9097832 

.9997707 

.0997733 

1  .9997633 

,  .9997482 

.9097251 

.0090867 

,9096102 

.9094208 

.9980362 

3.8 

.0098653 

.9098545 

.9908622 

.0998478 

.0998412 

1  .9998311 

.9998167 

>  .0997902 

.9997306 

.9990119 

.0992082 

3.0 

.0099038 

.9999033  ; 

.9999(118 

.0998989 

I  .0998945 

.9998878 

.0998776 

.9998608 

.9908276 

.9007420 

.9996020 

4.0 

.9990307 

.9099303 

.9990363 

.0999334 

.9999305 

.9999261 

.0990195 

.9909085 

.9908870 

.9998309 

.0000045 

4.1 

.9099587 

.9099585 

.9990578 

.9999500 

.9999547 

■  .0990610 

.0990475 

.9099404 

.9999200 

.9908900 

.9997763 

4.2 

.9999733 

.0099732 

.9099727 

.9990720 

.0999707 

.9990089 

.9009661 

.9009016 

.0999527 

.9999202 

.9908523 

4,3 

.0999829 

.9099828 

.9990820 

9090821 

.9099813 

9099301 

.9900783 

.0900754 

.9009008 

.9990548 

.0009034 

4.4 

.9009892 

.9999891 

.9999889 

.0990880 

.9990881 

.0999874 

.0090863 

.9099845 

.9090809 

.9099715 

.0000375 

4.5 

,0990032  i 

.9999932 

.9999931 

.9990929 

,0999926 

.9999921 

.9999914 

.9000002 

.0099881 

.0099822 

.0900600 

4.0 

.0999968 

.9999967 

.9999967 

.9990956 

.9999954 

.0999951 

.9090947 

.9900939 

,9999020 

.9999880 

.9900740 

4.7 

.9999074 

.9999974 

.9099973 

.9099973 

.0999971 

.9090970 

.9000967 

.9000003 

.9999055 

.0009932 

.0900840 

4.B  | 

.9999984 

.9999984 

.9999984 

.9099983 

.0990983 

.9990982 

.9999980 

.9990977 

.9999972 

.9999050 

.9909001 

4,1) 

.9999090 

.9999990 

.9999990 

.9990990 

.9999990 

.9999089 

.0099688 

.9990080 

.9099983 

,9000975 

.0990030 

6.0 

.0909004 

.9999994 

.9990094 

.9099904 

.9999994 

.9999993 

.9999993 

.0999092 

.9990990 

.9909085 

.0900063 

5.1 

.9909097 

.9999997 

.9999097 

.9099990 

.9090996 

.9999990 

.0699990 

.9900995 

.9999094 

.9999901 

.9990078 

6.2 

,0009998 

.9999998 

.9999998 

.9999998 

.0900098 

.9990098 

.9699998 

,9999997 

.9900097 

,0099905 

.9900087 

6.3 

. 9999690 

.9999999 

,9999999 

.9099999 

.9999099 

.9999990 

.9699090 

.0999998 

.9990098 

.9909997 

9999992 

3.4 

.0990099 

.9999999 

,9990999 

.9990991) 

.9099000 

.9999990 

.9999999 

.9999996 

.9999999 

,0900908 

.0099096 

6.6 

1.0000000 

1,0000000 

1,0000000 

1.0000000 

1.0000000 

1.0000000 

1.0000000 

.9999999 

.9999999 

.9990909 

,9909997 

6.6 

1 .0000060 

1.0000060 

.9999999 

.9999098 

6.7 

1.0000000 

.9990009 

6.8 

6.9 

6.0 

1.0000000 

PCX’.  t)*-Un  probability  that  ft  point  fftlU  luud*  ft  circle  *hoee  tutor  to  ftt  tl.o  ortoA*  ud  whoM  rftdiuft  to  K  \ lam  the  larger  etandanl  deviation,  c  being  the  ratio  cl  the  mailer 
•Uadftfd  deviation  to  the  larger  ttandard  deviation. 
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Table  111 

VAUIKS  OFK  C'.OHUKSl'OMDl N(J  TO  C.UMUI.AT'VK  PHOHAHI1.ITY  /' 


0,0 

0.1 

...... 

0.2 

0,3 

0.4 

0.5 

0.6 

0.7 

0.8 

0,9 

1.0 

.6000 

0.74993 

0.93365 

0.90621 

1 .05769 

1.11807 

1.17741 

.7600 

1 . 18473 

liil 

1.16246 

M4WU 

1.28534 

1 .35143 

1.42471 

1.60231 

1.58271 

1.66511 

.0000 

1 .94485 

1.04791 

1.65731 

1 .57383 

1.69918 

1.73708 

1 ,70152 

1.86253 

1.94761 

2.01236 

2.14597 

1.06253 

1 ,07041 

1 .08420 

2.03686 

2.14598 

2.23029 

2.33180 

2.44775 

.9750 

2.24140 

2.24365 

2.25053 

2.26255 

2.28073 

2.34581 

2.40356 

2.48494 

2.58090 

2.71620 

2.57689 

2.57778 

2 . 58377 

2.59421 

2.63257 

2.71515 

2.79009 

2.89743 

3.03485 

.0960 

2,80883 

2.81432 

2,83830 

2.85894 

2.88859 

2.03347 

3.00431 

3.11073 

3.25526 

.9976 

2.02500 

3.03010 

nrrf  i 

3.07144 

3.09871 

3.13969 

3.20580 

3.31099 

3.46164 

.9990 

3.29053 

3.29205 

3.29073 

3.33464 

3.35040 

3.36647 

3.45698 

3.55039 

3.71692 

Given: 

ol  *  15  meters 

c2  -  20  meters 

a  =  50 

For  the  computation  the  following  numbers  are  needed 


erf  **  225 


c r\  =  400 

sin 2  a  *  0. 5868 

Substituting  in  the  formula  for  cr * 


1 


2  x  0.5868 


[225  +  400  +  425 2  “  4  x  0.5868  x  225  *  400] 


I 

1.1736 


[625  +  v-390, 625  "  211,  248] 


0.852071625  h  ^179,377J 

0 . 852  [625  +  423] 

0.852  x  1048 
893 
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crt  -  1/093  e  29.9  me  ter  a 

o*  0.852  [625  “  423] 

-  0.852  x  202 

«  172 

note  tluit  the  numbers  are  the  same  as  lor  the  a2  eu 1 cu 1  at i on  except  for 
tbe  minus  sign. 

°y  "  *'172  -  13.1  meters 

t:  =  cry/at  »  13.1/29.9  =  0.438 

As  in  the  example  under  Method  1,  it  is  desired  to  determine  the  proba¬ 
bility  of  locution  of  tliia  point  within  u  circle  of  30  meters  radius, 

liurter* s  dutu  is  presented  in  terms  of  the  ratio  c ,  a  function  Kt  and 
probability.  The  function  K,  multiplied  by  the  larger  of  the  two  standard 
deviations  obtained  by  this  truns format  ion  method,  gives  the  value  of  the 
rudius  of  the  circle  of  the  corresponding  vulue  of  probability  shown  in  the 
table.  In  the  exumplc  here,  the  value  of  the  rudius  of  the  circle  for 

which  u  probability  is  desired  is  given  us  30  meters.  Thus  we  may  solve 

for  the  proper  value  ol  K  by  the  equation: 

Hud  i us  of  circle  *=  K 

30  *  K29.9 

K  «  1.003 

Then  on  graph  of  Tig.  11  (or  in  Table  11)  for  K  ~  1,0  and  e  =  0.44 

(interpolating)  read  P  =  0.62  ns  was  obtained  from  Method  1. 

An  alternative  presentation  of  Harter’s  data  is  given  in  Fig.  12 
and  Table  III  where  the  parameters  are  selected  so  as  to  provide  ready 
information  about  the  sizes  of  circles  of  specific,  probability  value 
associated  with  ellipses  of  varying  eccentricities.  These  are  convenient 
as  one  often  is  specifically  concerned  with  the  CEP ,  the  50%  probability 
circle  or  the  90%  circle. 
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FIG.  12  CONVERSION  OF  PROBABILITY  ELLIPSE  TO  CIRCLE 


2.3  GKOMKTRICAI.  KRROR  CONS  I  DKHATIONS 


From  the  information  that  can  be  derived  by  using  these  two  methods 
of  transformation  of  elliptical  error  duta,  one  may  develop  curves  which 
show  for  constant  values  of  initial  error  that  the  size  of  a  circle  of  a 
fixed  value  of  probability  varies  aa  a  function  of  the  angle  of  inter¬ 
section  of  the  lines  of  position. 

To  simplify  the  investigation  of  geometrical  factors,  it  is  initially 
desirable  to  consider  the  special  case  of  cTj  =  a,t  -  a.  Under  this  special 
condition,  the  long  formulas  for  cr t  and  may  be  drastically  simplified 
to  facilitate  computation  as  shown  below 


/2 


2  sin  ~  a 
0 


(<T,  =  cr2  ) 


/  2 


cr  « 
1 


2  cos  —  a 

2 


(cr. 


V 


Taking  the  ratio  of  these  two  values,  a  simple  formula 
ratio  c. 

cr  1 

y  1 

c  *  —  *  tan  -  a 

cr.  2 


is  found  for  the 


(Detailed  derivations  for  these  simplifications  may  be  found  in  Appendix  D.  ) 

Utilizing  these  simplified  formulas,  significant  parameters  of  error 
ellipses  have  been  tabulated  in  Table  IV  as  a  function  of  the  intersection 
ungle  a.  Using  the  CF.P  curve  of  Fig.  12,  values  of  the  CEP  have  been  cal¬ 
culated  for  each  angle,  showing  that  the  CEP  increases  as  the  angle  of 
intersection  decreases,  (The  tabulation  has  been  carried  out  only  for 
values  of  angles  less  than  90°— the  numerical  values  are  symmetrical  about 
this  value  of  angle.)  The  last  column  in  the  table  gives  the  factor  by 
which  the  CEP  for  angles  less  than  90°  is  greater  than  the  CEP  for  a  right 
utiglc.  This  magnification  of  error  curve  is  plotted  in  Fig.  13.  A  similar 
computation  has  been  performed  for  the  90%  probability  circle  as  it  may 
be  seen  that  the  curve  for  this  value  of  probability  has  a  slightly 
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ERROR  FACTOR 


Tuli I  o  IV 


SltiNl K  It '.ANT  I’AHAMKTKHS  OK 
KIIHOII  l.l.l. !  I'.SKS  WHEN 

°  1  "  °"2  "  1  ' 11 


a 

i 7 

X 

o 

> 

r 

n 

CUP 

KHHOH 

FACTOR 

<H) 

1.0 

1.0 

l  a  0 

1.177 

1.177 

1 . 00 

HO 

1.10 

0.024 

0 . 830 

1.078 

1.186 

1.01 

70 

i.2:i4 

O.HhS 

0.700 

0,  900 

1.228 

1.042 

(.0 

1.414 

0.1117 

0.  577 

0.914 

1 . 292 

1.099 

SO 

1.072 

0.782 

0.  466 

0.847 

1.420 

1.206 

45 

1.847 

0. 7o6 

0.414 

1 .  508 

1.281 

40 

2.00 

0.75:'. 

0.304 

0.7B3 

1.620 

1.076 

:io 

2.74 

0.733 

0.  208 

0.734 

2.  01 

1.710 

20 

4.0o 

0.71H 

0. 17(i 

0.700 

2.83 

2.  42 

10 

8.11 

0.710 

0.087 

0.68(1 

5.52 

4.69 

Krror  Factor  =  CKI1/ 1 . 1 77 


FIG.  13  CEP  MAGNIFICATION  v*.  INTERSECTION  ANGLE 
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Table  V 

90%  CIRCLE  ERROR  KACTOH 


differing  shape  from  the  CEP  curve — see 
Fig.  12.  Values  are  given  in  Table  V.  It 


a 

c 

K 

9  OX 

R 

ERROH 

FACTOR 

90 

1.0 

2. 145 

2.  145 

1.00 

80 

0.839 

1.98 

2.  18 

1.015 

70 

0.700 

1.86 

2.30 

1.07 

(t0 

0.  577 

1.775 

2.51 

1. 17 

SO 

0. 

1.72 

2.88 

1.34 

13 

0.414 

1.702 

3.15 

1.47 

40 

0.364 

1,687 

3.47 

1.615 

30 

0.268 

1.665 

4.53 

L.  11 

20 

0.  176 

1.652 

6,72 

3.13 

10 

0.087 

1.645 

13.35 

6.22 

Error  Faotiir  “  90S  fl/2.145 


is  well  known  in  any  problem  involving 
position  that  the  best  results  are  obtained 
when  the  crossing  angle  is  close  to  90°. 

The  curves  of  Fig.  13  indicate  the  magnitude 
of  the  growth  of  error  as  the  angle  varies 
from  90° 

it  is  also  of  interest  to  consider  an 
inverse  problem — what  values  of  probability 
result  if  the  radius  of  the  circle  iB  held 
constant  at  the  minimum  value  corresponding 
to  that  obtaining  for  the  intersection 


angle  =  90°?  An  answer  to  this  question  may  be  obtained  from  the  prob¬ 
ability  v  s .  intersection  angle  curves  given  under  Method  1,  Fig.  7. 


Along  the  ordinate  fi/cr  =  1.177  which  corresponds  to  the  CEP  for  the 
circular  case,  one  may  read  the  lesser  values  of  probability  corresponding 
to  the  various  intersection  angles.  Likewise 
one  may  also  obtain  the  probability  values 
corresponding  to  holding  a  circle  the  size  of 


Table  VI 


the  90%  probability  circle  for  the  circular. 


PROBABILITY  DECREASE  WITH 
ANGLE  FOR  A  CIRCLE  OF 


CONSTANT  RADIUS  o',  '  o2 


case  by  using  the  ordinate  fl/cr  =  2.15  (also 
equivalent  to  1,83  times  the  CEP).  These  two 
curves  are  plotted  in  Fig.  14  and  the  numer¬ 
ical  values  are  given  in  Table  VI.  It  is  to  be 


A  *  CEP/9Q0 
H  ~  1.83CKP/9Q0 
e  90*  Probability  at  90° 


noted  that  the  probability  values  ure  not  in¬ 
versely  related  to  the  error  factors  plotted 
in  the  preceding  curves.  The  geometric  error 
factor  was  shown  to  be  a  simple  trigonometric 
function;  the  probability  curves  are  exponen¬ 


a 

P 

P 

90 

50 

90 

80 

49.4 

89.2 

70 

47.5 

86.9 

60 

44,0 

82.4 

30 

39.5 

76 

40 

37 

66 

30 

25 

53 

20 

17 

37 

10 

8 

19 

tial  functions. 
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FIG.  14  DECREASE  IN  PROBABILITY  FOR  A  CIRCLE  OF  CONSTANT 
RADIUS  v*.  ANGLE 
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3.  ANALYSIS  OF  MULTIPLE  ERROR  ELLIPSES 


The  tools  developed  in  the  preceding  section  permit  the  specifica¬ 
tion  of  individual  error  ellipses  about  a  single  point.  As  shown  in 
Fig,  1  a  rcul  problem  in  pos i t ion- 1  oca t ion  involves  the  consideration  of 
the  combination  of  errors  from  a  number  of  sources.  And  us  the  preceding 
section  showed,  in  generul  each  of  these  various  sources  of  error  will  be 
expressed  as  an  error  ellipse.  Following  the  methods  of  the  previous 
section,  each  ellipse  can  be  expressed  in  terms  of  the  standurd  deviations 
along  its  major  and  minor  axes,  The  problem  of  the  combination  of  multiple 
error  ellipses  is  the  determination  of  the  proper  method  combining  a 
number  of  individual  errors  to  obtain  the  total  error  at  some  desired 
point.  In  the  general  case,  the  ellipses  will  not  be  oriented  relative 
to  one  another  in  any  wuy  hut  a  random  manner.  The  exact  consideration 
of  this  random  orientation  of  the  axes  of  the  ellipses  complicates  the 
analysis,  but  it  is  necessary  to  ubtain  accurate  answers.  More  facile, 
but  approximate  methods,  sometimes  seen  in  the  literuture,  will  be  dis¬ 
cussed  at  the  end  of  this  section.  Such  approximate  methods  often  result 
in  sizuble  errors — errors  which  almost  always  come  oiit  in  the  wrong 
direction  so  that  the  system  appears  to  be  better  than  it  really  is. 

3  .1  SPECIAL  CASE  — MUTUALLY  PARALLEL  AXES 

Before  attempting  the  analysis  of  the  generul  cuse,  it  is  helpful 
to  look  at  the  restricted  und  unlikely  case  where  all  the  various  ellipses 
to  be  considered  have  their  uxes  mutuuily  parallel.  This  special  case 
will  then  leud  to  the  more  generul  cuse. 

Referring  again  to  Fig.  I,  there  urn  four  error  ellipses  of  interest: 

Weapon  dispersion 

Gun  location 

Forward  observer  location 

Target  location  with  respect  to  forward  observer. 

In  a  system  with  these  errors  we  wish  to  obtain  the  probability  of  damage 
to  the  target  assuming  that  the  shell  must  land  within  a  circle  of 
20-meter  radius  in  order  to  obtain  the  desired  damage  level. 
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For  this  initial  exumple  it  is  assumed  that  ull  the  error  ellipses 
have  their  uxes  mutually  pn ra 1 1 e l — ul igncd  north  and  east,  for  example — 
and  have  the  parameters  listed  in  the  example  shown.  Standard  deviations 
for  the  two  uxes  of  each  ellipse  are  listed.  The  method  of  obtaining  the 
total  error  at  the  target  is  to  obtain  the  sum  of  the  variances  in  the 
two  directions  and  to  convert  these  sums  to  the  standard  deviations  of  the 
total  error  ellipse  ul  the  turget.  The  desired  probability  may  then  be 
obtained  by  Method  2  of  the  preceding  section.  The  calculations  are 
s  liown  in  do  tu  i  l . 


GIVEN 

STANDARD  DEVIATIONS 

cr  ov 

«  y 

VARIANCES 

*7“P“ 

*  y 

Weapon  dispersion 

3  m  40  m 

9 

1600 

Gun  location 

10  15 

100 

225 

Forward  observer 

15  20 

225 

400 

Target  location 

30  10 

900 

100 

Sums  of  variances 

1234 

2325 

Take  square  roots  of  variances  to  ubtain  new  standard  deviations  of  total 
error  ellipse  ut  target. 

»  /l234  »  35.1  m 
«  v/2325  •  48.2  m 

From  Method  2,  Section  2‘ 


c 

Hudius  of  circle 


20 

K 


35.1/48.2  ’  0.729 

Kcr. 

1  •  r  g  «  r 

48.2  K 

20/48.2  *  0.415 


Then  from  the  graph  of  Fig.  11  i'or  K  =  0.415  and  c  =  0.73  (interpo¬ 
lating)  read  P  «  0.11. 


In  this  case  of  the  combination  of  several  ellipses  cr  has  turned 
out  to  be  larger  than  cr  .  In  such  cases  the  factor  K  is  always  to  be 
multiplied  by  the  larger  of  the  two  values  of  sigma,  cr^  or  cr  to  obtain 
the  radius  of  the  probability  circle.  The  formulas  given  on  Fig.  12  for 
simplicity  are  stated  in  terms  of  cr  always  being  the  larger.  The 
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formulas  ussure  this  condition  for  uny  single  ellipse  taken  uloije  without 
reference  to  other  ellipses.  However,  when  ellipses  uru  combined,  either 
stundard  deviation  may  turn  out  to  be  the  lurger. 

By  way  of  comparison  it  is  interesting  to  calculate  how  much  of  the 
error  is  contributed  by  location  measurement  errors  and  how  much  is  con¬ 
tributed  by  the  dispersion  of  the  weapon.  If  we  consider  the  weapon 
dispersion  ellipse  alone 


c 

Radius  of  circle 


20 

K 


3/40  =  0.075 

Kcr. 

I*r|tr 

40  K 
0.5 


Then  on  the  graph  of  Tig.  11  for  K  =  0.5  and  c  *  0.075  (interpolating) 
read  P  e  0. 37. 

Thus  with  perfect  location  of  all  elements,  this  gun  with  the  stated 
dispersion  (105  mm,  mid-range)  has  a  37%  probability  of  landing  a  shell 
within  a  circle  of  20-meter  radius.  But,  when  the  three  location  errors 
are  combined  with  the  dispersion,  the  probability  falls  to  11%. 

The  numbers  shown  are  realistic  for  a  system  of  good  accuracy  and 
better  than  most  performance  of  today. 


This  method  of  adding  variances  along  the  two  axes  at  right  Bngles, 
down  range  and  cross  range,  is  the  standard  method  of  preparing  an  error 
budget  for  a  weapon  system.  The  method,  however,  is  not  sufficient  when 
one  wishes  to  combine  ellipses  having  random  orientations  of  their  axes. 
And  since  statistical  distributions  are  involved,  simple  trigonometric 
resolutions  from  one  set  of  axes  to  another  are  insufficient. 


3.2  GENERAL  CASE— RANDOM  ORIENTATION  OF  AXES 

In  the  general  case  of  random  orientation  of  axes  in  any  number  of 
error  ellipses,  a  more  complex  procedure  for  combination  is  necessary 
than  that  described  in  the  previous  subsection.  Briefly,  a  reference 
set  of  axes  must  be  chosen  and  the  orientation  of  each  error  ellipse 
with  respect  to  these  axes  must  be  determined.  Then  the  variances  along 
these  axes  must  be  computed,  a  procedure  which  will  also  involve  cross 
product  terms  (See  Appendix  B).  The  special  variances  and  the  cross 
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product  functions  muy  then  bo  added  to  obtain  tlio  correspond  i  ng  (unctions 
of  the  final  ellipse.  From  these  two  vuriunoos  und  tlio  cross  product 
function  of  the  final  ellipse,  which  are  ussociuLod  with  the  urhi Lrnri I y 
chosen  set  of  axes,  final  values  of  and  'J  ulong  the  major  and  minor 
axes  of  the  ellipse  may  be  calculated.  Necessary  formulas  are  givun 
below.  Illustrative  examples  will  clarify  the  description.  The  labor 
involved  is  considerable  for  any  real  example.  lor  this  reason,  a  com¬ 
puter  program  was  developed  to  permit  such  calculations  to  be  run  off  in 
largi  quantities  as  a  part  of  the  systems  evaluation  program.  Derivations 
of  the  formulas  presented  are  to  be  found  in  Appendix  D. 

An  example  is  illustrated  in  Fig.  15.  The  three  smaller  ellipses 
are  the  given  inputs  to  the  problem  and  the  large  ellipse  on  the  right 
represents  the  combination  of  the  three  smaller  ones.  Each  of  the  three 
small  ellipses  is  described  in  terms  of  its  own  cr^  and  cr^ .  Also  shown  are 
the  angles  between  the  x-axis  of  each  ellipse  and  the  arbitrarily  selected 
reference  axes  which  are  designated  the  w  and  i  axes.  To  obtain  the 
parameters  of  the  final  ellipse,  variances  for  each  ellipse  along  the 
id  and  z  axes  are  calculated.  Because  the  axes  are  not  those  of  the 
individual  ellipses,  an  additional  function  p-,  the  cross  product  function 
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FIG.  15  ELLIPSES  WITH  RANDOM  ORIENTATION  OF  AXES 
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is  also  required  for  each  ellipse.  These  three  functions,  the  two  vari¬ 
ances  and  the  cross  product  function,  from  each  ellipse  aro  then  added 
to  obtain  the  corresponding  functions  in  the  final  ellipse.  Functions 
lor  the  final  ellipse  are  designated  with  the  second  subscript  f.  These 
three  functions  in  the  fins  ellipse  may  then  be  converted  to  the  a  * j 
and  v ^  of  that  ellipse  along  its  major  and  minor  axes.  The  calculations 
are  tabulated  for  this  illustrative  example,  on  the  next  few  pages,  fol¬ 
lowing  the  listing  of  the  necessary  formulas.  The  problem  illustration 
gives  all  intermediate  values  needed  but  does  not  show  every  individual 
calculation  required.'  The  extensive  amount  of  work  necessary  for  such 
a  simple  problem  clearly  shows  the  need  for  a  mechanized  computer  solution. 
The  simple  problem  illustrated  requires  an  hour  or  more  of  hand  labor  with 
tables  and  slide  rule-— the  computer  handles  it  in  three  seconds. 

The  necessary  formulas  for  this  calculation  are  given  below  without 
derivations.  Further  discussion  and  derivations  will  be  found  in  Appendix  E. 
The  formulas  are  given  in  general  notation  using  the  letter  i  to  represent 
the  general  ellipse.  In  the  use  of  the  formulas,  i  =  l,2,3,.,,n  according 
to  the  number  of  ellipses  involved.  The  formulas  are  mo3t  conveniently 
expressed  in  terms  of  some  auxiliary  functions  which  are  computed  first. 


Define 


Then 
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We  now  have  the  parameters  of  the  final  ellipse  in  terms  of  the  w 

and  2  axes,  To  eliminate  the  cross  product  function  p-.  and  obtain  cr 

/  *  / 

cr  ^  along  the  major  and  minor  axes  of  the  final  ellipse,  we  use  the 
formulas  below. 


The  orientation  of  the  final  ellipse  with  respect  to  the  w  and  i 
axes  is  given  by  the  formula  below. 

tan  26  ^ 

The  numerical  example  follows. 


GIVEN 

CONDITIONS 

ELLIPSE  H 

ELLIPSE  *2 

ELLIPSE  #3 

a 

V!  1 

15  m 

10  m 

10  m 

10  ra 

20  m 

20  m 

0 

i 

45° 

60° 

150° 

calculated 

2 

axi 

225 

100 

100 

2 

ayi 

100 

400 

400 

36 


37 


crf  I  =  v' 600, 3  =  24.5  m 

„  2  x  0.09441  x  25.74  *  25.74 

ton  2  6.  3  - - - 

f  662.9  “  662.9 

numerator 

a  «-  .-1  -i 

0 

"  oo 

arcton  c a  5  90" 

2  6{  =  .  90° 

=  45° 

•  The  computer  solution  of  this  same  problem  gives  answers  differing 
only  by  one  in  the  fourth  place,  thus  confirming  the  results. 

Because  of  the  labor  involved  in  the  foregoing  calculations,  some 
references  have  indicated  that  errors  may  be  combined  by  converting  each 
error  distribution  to  the  CEP  at  that  point  according  to  Method  2  and 
then  combining  the  individual  CEP  values  root»sum-square.  This  method 
almost  always  has  been  found  to  give  too  small  a  value  of  the  CEP  at  the 
final  ellipse.  Conversion  to  individual  circles  eliminates  the  effect 
of  orientation  of  the  ellipse  which  is  an  important  consideration  in 
combination.  To  illustrate,  the  same  three  ellipses  used  in  the  pre¬ 
viously  lengthy  example  will  be  handled  by  this  simplified  method  and 
compared  with  the  answer  obtained  by  converting  the  final  ellipse  to  a 


ELLIPSE  *1 

ELLIPSE  *2 

ELLIPSE  #3 

Or 

X  i 

15 

10 

10 

O  . 

7* 

10 

20 

20 

c 

0.667 

0.50 

0.50 

K 

0.956 

0.870 

0.870 

CEP 

14.35 

17.4 

17.4 

CFP2 

206 

302 

302 

Sum  of 

(CEP1 a) 2 

206 

302 

302 

810 


CEP  of  final  location  =  810  “  28.45  m 
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Now  compare  this  with  CEP  of  final  ellipse  as  just  computed 


cr  ,  =  26.93 

*/ 

<r  ,  ■  24.5 

7  / 

c  -  24.5/26.93  =  0.948 

K  -  1.144 

CEP  -  L. 144  x  26.93  =  30.8  m 

The  error  in  this  case  is  relatively  small  while  the  simplified 
method  gives  an  answer  that  is  8%  too  small.  Other  sample  calculations 
indicate  that  errors  of  up  to  20%  can  occur  using  this  simplified  method 
Approximate  calculations  of  the  CEP  are  discussed  in  more  detail  in 
Appendix  A. 
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APPENDIX  A 

MEASURES  OF  ERROR 


41 


APPENDIX  A 


MEASURES  OF  ERROR 


1.  INTRODUCTION 

This  appendix  describes  the  various  terms  used  as  measures  of  error 
in  more  detail  than  the  discussion  in  the  body  of  this  Memorandum.  In 
addition,  other  terms  often  noted  in  the  pertinent  literature  are  defined 
and  reluted  to  the  terms  used  in  this  report. 

The  reuder  is  cautioned  to  reud  analyses  of  systems  accuracy  with 
care.  The  literature  examined  by  the  SRI  project  team  in  the  progress  of 
this  study  sometimes  contains  numerical  errors.  More  serious,  however, 
incomplete  or  misleading  definitions  of  terms  ure  often  found,  and,  occa¬ 
sionally,  incorrect  definitions  have  been  noted.  In  other  references  it 
may  be  impossible  to  tell  exuctly  what  whs  meant  by  a  particular  measure 
of  error.  The  discussion  of  measures  of  error  in  this  appendix  is  in¬ 
tended  to  help  clour  uwuy  such  misunderstanding  and  confusion. 

2.  ONE  DIMENSIONAL  ERROR  TERMS 

Although  the  basic  problem  of  position  locution  is  concerned  with  the 
two  dimensions  necessary  to  describe  an  ureu,  one  dimensional  error  mea¬ 
sures  ure  commonly  applied  to  each  of  the  two' dimensions  involved.  In 
fact,  as  shown  in  the  discussion  in  the  body  of  this  memorandum,  it  is 
most  convenient  to  do  this  to  permit  u  truly  general,  approach  to  the  con¬ 
sideration  of  error  ellipses.  Thus  the  meusures  of  error  concerned  with 
one  dimensional  of  Gaussian  distributions  are  important.  The  following 
terms  ure  frequently  met  and  each  is  described  in  following  paragraphs: 
standard  deviation,  RMS  error,  signm,  probable  error,  and  variance. 

2.1  Stan  da  Hi)  Deviation,  RMS  Eh  non,  Siuma  io-> 

The  three  terms  all  mean  the  same.  The  basic  equation  of  the  normal 
distribution  indicates  the  use  of  the  Greek  letter  sigma  from  which  the 
shorthand  use  of  sigma  foi  standard  deviation  arises. 


43 


Preceding  page  blank 


Standard  deviution  of  a  measurement  system  is  a  property  that  may  be  de¬ 
termined  experimentally.  If  u  large  number  of  measurements  of  the  same 
quantity— *u  length,  for  example  — are  mode  und  compared  with  a  standard, 
the  standard  deviution  is  the  square  root  of  the  sum  of  the  squares  of 
the  deviations  from  the  meun  or  uveruge  value  divided  by  the  number  of 
measurements  taken.  Symbolically,  this  operation  is  represented  as 


o 

The  term  rms  ( root -mean -square  )  error  comes  from  this  latter  method. of 
computation. 

Numerically,  one  sigmu  corresponds  to  68%  of  the  distribution!  that 
is,  if  u  large  number  of  measurements  were  made  of  a  given  quantity,  68% 
of  the  errors  would  be  no  gt'ea to r .  than  the  valun  of  one  standard  deviation 
Likewise  2 o  corresponds  to  95%  of  the  total  errors,  and  3cr  to  99.6%  of 
the  total  errors. 

2.2  Phobabi.K  Ehhor 

This  term  is  identical  in  concept  to  standard  deviation  when  con¬ 
sidered  as  the  rms  error  determined  after  a  series  of  measurements.  The 
term  differs  from  standard  deviation  in  that  it  refers  to  that  value 
corresponding  to  the  median  error;  no  more  titan  half  the  errors  in  the 
measurement  sample  are  greater  than  the  value  of  the  probable  error. 

Linear  probable  error  is  related  to  standard  deviation  by  a  multiplying 
factor,  One  probable  error  equals  0,6745  times  one  standard  deviation. 
Probable  error  is  used  in  Army  artillery  manuals  as  a  measure  of  weapon 
component  errors,  such  as  the  range  and  deflection  errors  associated 
with  u  particular  weapon  and  ammunition.  Industrial  practice  in  the 
United  Stales  also  employs  the  probable  error  when  a  measurement  is  re¬ 
reported  in  the  manner  of  173.23  ±  0.05  ft.  The  0,05  ft  is  to  be 
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interpreted  as  a  probable  error  of  the  measurement.  Probable  error  has 
not  been  used  in  the  analyses  of  this  study  for  it  is  not  as  convenient 
un  error  measure  to  handle  mathematically  as  the  standard  deviation. 

2.3  Variance 

This  term  is  met  most  frequently  in  detailed  mathematical  discus¬ 
sions.  The  term  refers  to  the  square  of  a  standard  deviation.  It  is 
useful  in  simplifying  the  algebra  of  some  complex  mathematical  deriva¬ 
tions  (see  Appendix  B  for  examples).  It  also  is  a  convenient  concept 
when  preparing  an  error  budget  made  up  of  many  separate  components  of 
error  for  the  individual  variances  may  be  added  directly  to  obtain  the 
total  variance. 

3.  TWO-DIMENSIONAL  ERROR  TERMS 

Terms  similar  or  identical  in  words  to  those  used  for  one-dimensional 
error  descriptions  are  also  used  with  two-dimensional  or  bivariate  error 
descriptions.  However,  in  the  two-dimensional  case  ,  not  a  1 1  of  these  terms 
have  the  same  meaning  as  before  and  considerable  care  is  needed  to  avoid 
confusion. 

3.1  Standard  Deviation  or  Sioma 

These  two  terms,  used  interchangeably,  have  a  definable  meaning  only 
ill  the  specific  case  of  the  circular  normal  distribution  where  <7  x  =  cr  y 


In  the  case  of  the  circular  normal  distribution,  the  standard  deviation 
a  is  equivalent  to  the  standard  deviation  along  both  orthogonal  axes. 
Because  we  are  here  concerned  with  a  radial  distribution,  the  total  dis¬ 
tribution  of  errors  involves  different  numbers  from  those  of  the  linear 
case.  In  the  circular  case,  1  a  error  indicates  that  39.3%  of  the  errors 
would  not  exceed  the  value  of  the  1  c r  error;  86.4%  would  not  exceed  the 
2  o-  error,  and  98.8%  would  not  exceed  the  3  cr  error. 

Because  the  usual  case  where  there  are  two-dimensional  distributions 
is  that  the  standard  deviations  along  the  two  axes  are  different,  resulting 
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ill  mi  elliptical  distribution,  the  circular  standurd  deviation  is  less 
useful  thun  the  linear  stundard  deviation.  It  is  more  common  to  describe 
two-dimensional  distributions  by  the  two  separate  one -dimens i ona l  stan¬ 
dard  deviations  associated  with  ouch  error  axis,  Ileferences  often  do  not 
make  this  distinction,  however,  referring  to  the  position  accuracy  of  u 
system  as  "pOO  ft  (2  :>■),"  for  example.  Such  a  description  leaves  the 
reader  wondering  whether  the  measure  is  circular  error,  in  which  case 
the  numbers  describe  the  86%  probability  circle,  or  whether  the  numbers 
are  to  be  interpreted  ns  one  -d i mens i ona 1  signals  along  each  axis,  in  which 
case  the  95%  probability  circle  is  indicated  (assuming  the  distribution 
to  be  circular,  which  actually  it  may  not  be).  The  analyses  of  this  re¬ 
port  have,  in  general,  used  the  two  separate  linear  standard  deviations 
as  error  measures.  Where  specific  circular  measures  have  been  used,  they 
are  so  noted  carefully  to  avoid  confusion.  (See  next  subsection.) 

The  term  RMS  error  when  upplied  to  two rd imens i ona 1  errors,  does  not 
have  the  same  meaning  as  the  standard  deviation.  The  term  is  often  used 
in  the  literature,  although  it  has  an  ambiguous  meaning  in  relation  to 
terms  of  probability.  For  this  reason  its  use  hus  been  deprecated  in 
this  study.  The  term  is  discussed  separately  in  a  following  subsection. 

3.2  CEP  (Cihculan  Error  Probable —Also  Sometimes 
CPE,  Circular  Probable  Ehhoh) 

lit  a  circular  normal  distribution  this  term  refers  to  the  radius  of 
a  circle  containing  50%  of  the  sample  of  the  individual  measurements  be¬ 
ing  made,  or  the  radius  of  the  circle  inside  of  which  there  is  a  50% 
probability  of  being  located.  This  Is  a  common  measure  often  used  with 
weapon  systems  and  position  location  systems. 

The  term  CEP  is  also  used  to  indicate  the  radius  of  a  circle  inside 
of  which  there  is  a  50%  chance  of  being  located,  even  though  the  actual 
error  figure  is  an  ellipse  (Fig.  A-l).  The  body  of  the  Memorandum  de¬ 
scribes  the  method  of  obtaining  such  CEP  equivalents  when  given  ellipses 
of  varying  eccentricities.  Curves  and  tables  are  furnished  to  perform 
this  calculation.  In  the  literature,  despite  the  availability  of  these 
curves  and  tables,  approximations  are  often  made  for  this  calculation  of 
a  CEP  when  the  actual  error  distribution  is  elliptical.  Several  of  these 
approximations  are  indicated  and  plotted  for  comparison  with  the  exact 
curve  in  Fig.  A-2.  Of  the  various  approximations  shown,  the  top  curve, 
the  one  which  diverges  the  most  rapidly,  appears  to  be  the  most  commonly 
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FIG.  AO  ERROR  ELLIPSE  AND  CIRCLE  OF  EQUIVALENT  PROBABILITY 

used  in  published  systems  unulyaes.  Use  of  the  curves  and  tables  given 
in  the  body  of  this  Memorandum  is  recommended  to  avoid  such  problems  of 
approximations, 

Another  factor  of  interest  concerning  the  relationship  of  the  CEP  to 
various  ellipses  is  that  the  areu  of  the  CEP  circle  is  always  greater  than 
the  basic  ellipse.  Calculations  made  using  the  values  of  the  tables  given 

in  the  body  of  the  Memorandum  are 
given  in  Table  A-l  where  it  may  readily 
be  seen  that  the  divergence  between  the 
uctual  urea  of  the  ellipse  of  interest 
and  the  circle  of  equivalent  probability 
increases  os  the  ellipse  becomes  thin* 
ner  and  more  elongated.  This  fact 
provides  a  powerful  reason  for  the 
method  of  analysis  used  in  this  study 
of  considering  the  ellipses  directly, 
especially  when  a  number  of  ellipses 
are  involved  in  the  determination  of 
a  final  probability  figure. 


Table  A-l 

COMPARISON  OF  AREAS  OF  50%  ELLIPSES 
OF  VARYING  ECCENTRICITIES  WITH  AREAS 
OF  CIRCLES  OF  EQUIVALENT  PROBABILITY 


C  ■  a/t 

AREA  OF 

50%  ELLIPSE 

AREA  OF 

EQUIVALENT  CIRCLE 

0.0 

0 

1.43 

0.1 

0,437 

1.46 

0.2 

0.874 

1 . 56 

0.3 

1.31 

1.76 

0.4 

1.75 

2.06 

0,5 

2.08 

2.37 

0.6 

2.62 

2.74 

0.7 

3.05 

3.12 

0.8 

3.49 

3.52 

0.9 

3.93 

3.94 

1.0 

4,37 

4.37 
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FIG.  A-2  CEP  FOR  ELLIPTICAL  ERROR  DISTRIBUTION  APPROXIMATIONS 
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Tin*  effect.  of  con  vers  i  on  of  ellipses  to  t,  h  «•  equivalent  CEP  values 
is  well  illustrated  in  t.lic  nimbi  nut  ion  of  ellipses  of  different  orictitu- 
liniis.  If  two  equal  ellipses  whose  ruin  live  or  i  ontut  i  on  in  90"  are  com¬ 
bi  noil,  the  result  is  u  circle.  Ah  mentioned  in  the  lust  section  of  the 
body  of  the  Memorandum,  combination  of  ellipses  is  sometimes  calculated 
by  converting  eueli  individually  to  its  equivalent  CKP,  and  then  combin¬ 
ing  the  individual  CKP  values  rim  t. -sum -squn  re .  The  simple  example  shown 
hero  indicates  that  a  considerable  degree  of  error  cun  result  from  such 
a  nimbi  nut  ion  of  separately  obtained  CKP  values  when  compared  with  the 
CKP  obtained  .from  the  combined  figure. 


ELLIPSE  m  ELLIPSE  ff 2 


O'  ,  ■ 

1 

<7  = 

L0 

»  1 

*  2 

O'  ,  = 

y  1 

10 

cr  .  ■ 

y  l 

1 

C  = 

0.1 

C  = 

0.1 

K  - 

0.681 

K  = 

0.681 

CKP  - 

Kcr,i 

CEP  - 

In  each  case  CEP 

=  6.81.  If 

the  two 

values  of  CEP  are 

combined  root -sum 

-square  the 

result 

is  6,81  /2*  *  9.65. 

If  we  now  combine 

variances 

to  obtui 

n  the  combination 

the  results  are: 

cr2  .  « 

1 

a2,  ’ 

100 

i  1 

y  l 

*lt  - 

100 

1 

101 

cr2  .  « 

y  f 

101 

a«/  = 

10.04 

ayf  “ 

10.04 

Thus,  the  final 

figure  is  a 

circle. 

C 

=  1.0 

K 

-  1.177 

CEP 

=  1.177  * 

10.04 

-  11.9. 
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For  the  75%  probability  circle 


1 


2<r‘ 


0.75 


R2 

2a2 


R  - 

R  = 

fl7  5 
fl50 


0.25 


In  0.25 


In  4 


cr  /2  in  4 

1 . 66  5o" 

1.66.5u- 
1 , 1 7  7i/ 


1.414 


The  factors  tabulated  in  Tuble  A-2  are  sometimes  used  in  the  litera¬ 
ture  to  relate  varying  probability  circles  when  the  basic  distribution  is 
not  circular,  but  elliptical.  That  such  a  procedure  is  inuccurote  maybe 
seen  by  the  curves  of  Fig.  A-4.  These  curves  were  prepared  from  the  val¬ 
ues  of  Table  III  in  the  body  of  the  Memorandum.  It  may  be  seen  that  the 
errors  involved  are  smull  when  small  e 1 1 i pti ci t i es  ure  involved.  But  the 
errors  increase  significantly  when  both  high  vulues  of  probability  are 
desired  and  when  the  ellipticity  increases  in  the  direction  of  long,  nar¬ 
row  distributions. 


3.3  Comu;  lai  ion  Coki- r  i c ikn t 

In  many  statistical  references  the  presence  of  u  cross-produc t  term 
is  indicated  as  a  correlation  between  the  two  variables,  such  as  z^andzj 
As  used  in  this  statistical  sense  the  term  correlation  does  not  imply 
that  the  variables  Zj  and  n 2  are  connected  in  any  physical  way.  The  two 
variables  are  independent  and  a  change  in  one  does  not  affect  the  other. 
The  factor  p  (rho)  which  appears  in  equations  relating  multiple  ellipses 
to  an  arbitrary  set  of  axes  is  that  factor  often  called  the  correlation 
coefficient  in  statistical  literature.  Because  correlation  cun  be 
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FIG.A-4  RELATION  OF  PROBABILITY  CIRCLES 
TO  CEP  vs.  ELLIPTICITY 


misinterpreted  ii.H  implying  u  physical  i  liter  re  lu  t  i  oush  i  p  wliich  does  not. 
actually  exist,  the  use  of  the  term  bus  linen  avoided  in  this  study,  de¬ 
spite  the  consequent  need  for  the  circumlocution  of  11  cross-product  term 
in  z  ,z,.  " 

The  term  correlation  is  also  used  in  some  of  the  references  to  indi¬ 
cate  the  case  when  there  is  an  interrelationship  between  the  two  variables, 
a  condition  that  ohtuins  in  some  physical  systems  where  two  signals  may 
be  synchronized  to  a  third.  Hence,  one  must  he  very  careful  when  using  the 
term  to  insure  n  proper  understanding  of  the  mathematics.  Because  the 
term  correlation  is  thus  subject  to  possible  confusion  and  con  i.  rad  i  c  t  i  o.-  , 
its  use  when  referring  to  a  cross-product  term  has  been  avoided  in  this 
study.  When  it  occurs,  it  will  refer  solely  to  an  uutuui  physical  inter¬ 
relationship. 

3.4  Radial  oh  RMS  Khhoh,  d 

r  m  s 

The  terms  radial  error,  RMS  error,  and  d  are  identical  in  meaning 
when  applied  to  two -ditnensi oau  1  errors.  Figure  A-5  illustrates  the  defi¬ 
nition  of  d  ,  Tt  is  seen  to  be  the  square  root  of  the  sum  of  the  squures 
of  the  one  sigma  error  components  along  the  mujor  and  minor  axes  of  a 
probability  ellipse.  The  figure  details  the  definition  of  l  d  .  Sim¬ 
ilarly,  other  values  of  d  may  he  derived  by  using  the  corresponding 
vulues  of  sigtnu.  The  measure  d(  is  nvt  equal  to  the  squure  root  of 
the  sum  of  the  squures  of  the  cr(  uud  c2  that  ure  the  busic  errors  con¬ 
nected  with  the  lines  of  position  of  a  particular  position  locution  sys¬ 
tem.  The  procedures  described  in  the  main  portion  of  this  Memorandum 
and  derived  in  Appendix  D  must  first  he  utilized  to  obtain  the  values 
shown  as  and  c r^  in  Fig.  A-5.  The  three  terms  used  as  a  measure  of 
error,  RMS  error,  rudiul  error,  and  d  ure  .somewhat  confusing  because 
they  do  not  correspond  to  u  fixed  vuluc  of  probability  for  a  given  value 
of  the  error  measure.  The  terms  cun  be  conveniently  re  luted  to  other 
error  measures  only  when  c t  •  cr  und  the  probability  figure  is  a  circle. 

In  the  more  common  elliptical  cases,  the  probability  associated  with  a 
fixed  value  of  d(  varies  as  a  function  of  the  eccentricity  of  the 
ellipse,  1  df  is  defined  us  the  radius  of  the  circle  obtained  when 
-r ,  “  1  in  Fig.  A-5  and  cr  vuries  from  zero  to  one.  Likewise,  2  d  i 
the  radius  of  the  circle  obtuined  when  cr x  -  2  and  tr  vuries  from  zero 
to  two.  Values  of  the  length  of  the  radius  d  may  be  calculated  us 
shown  mi  Table  A-3.  From  these  values  the  ussociuted  probabilities  may 
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FIG.  A-5  ILLUSTRATION  OF  ROOT  MEAN  SQUARE  ERROR 


be  dutermi ned  from  the  tubles  of  tlie 
body  of  this  annex.  The  variations 
of  prohubi lity  associated  with  the 
vulues  1  J  und  2  d  are  shown 

r  »n  s  r  in  8 

in  the  curves  of  Figs,  A-6  u mi  A-T, 

Fig.  A-8  shows  the  luck  of  u  constant 
relationship  in  u  slightly  different 
way.  Here  the  ratio  d  /CEP  is 
plotted  against  the  sume  measure  of 
ellipticity.  The  three  figures  show 
graphically  that  there  is  not  a  con- 
stunt  vulue  of  probability  associated 
with  a  single  vulue  of  d  .  While 
this  variation  is  not  great,  rt  is 
felt  to  be  unneeessuri  ly  confusing.  Thus  the  measure  d  Inis  not  been 
used  in  this  study. 
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cr 

y 

cr 

X 

Lbi NOTH  OK 

1  d 

r  ms 

PROBABILITY 

1  d 

v  ms 

2  d 

r  h.s 

0.0 

1.0 

1 . 000 

0.6U3 

0.954 

0. 1 

).() 

1 ,  (IDS 

0.682 

0.955 

0.2 

1.0 

1.020 

0.6H2 

0.957 

0.3 

1.0 

1 . 042 

0.676 

0.961 

0.4 

1.0 

1.077 

0.671 

0.966 

0.3 

1.0 

1.  1  111 

0.662 

0.9(i  9 

0.6 

1.0 

1 ,  1 66 

0.650 

0.973 

0.7 

1.0 

1 . 220 

0.641 

0.977 

O.H 

1.0 

1.280 

0.635 

0.9110 

0.9 

1.0 

1 . 345 

0.632 

0.9H1 

1.0 

1.0 

1.414 

0.632 

0.982 

/  ■’> - r* 

d  a  Ycr‘  +  a  winner  and  cr  writ  «t  right. 
hub  x  y  x  y 

u ii ^  Ian  tu  euuli  oilier* 
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APPENDIX  B 


DERIVATION  OF  METHOD  2  FORMULAS 

1.  INTRODUCTION 

In  the  body  of  this  memorandum  it  wus  shown  that  the  most  general 
case  of  analyzing  the  accuracy  of  a  positioning  system  resulted  in  the 
consideration  of  basic  input  error  vectors  intersecting  at  any  angle. 
Graphs,  nomograms,  and  tables  were  shown  that  may  be  used  to  obtain 
solutions  to  problems  of  accuracy  determination  where  the  lines  of  posi¬ 
tion  do  not  intersect  at  right  angles.  Two  methods  of  solutions  were 
given  there.  This  appendix  presents  a  complete  derivation  of  the  for¬ 
mulas  to  obtain  the  standard  deviations  along  the  major  and  minor  axes 
of  un  error  ellipse,  These  functions  are  then  used  as  the  inputs  to 
Harter’s  method  of  determination  of  elliptical  probabilities  described 
at  the  end  of  the  Method  2  discussion. 

The  presentation  of  the  Method  2  analysis  is  given  ahead  of  that 
for  Method  1  (see  Appendix  D)  because  of  the  detailed  discussion  developed 
during  this  analysis,  which  is  basic  in  philosophy  to  both  methods.  Also, 
some  of  the  formulas  developed  in  this  appendix  arc  useful  in  derivations 
of  some  of  the  intermediate  Method  1  solutions. 

Nonorthogona 1  bivariate  distributions  arise  when  the  lines  of  posi¬ 
tion  do  not  intersect  at  right  angles.  Such  a  condition  is  the  usual 
condition  whatever  the  type  of  navigational  system  used — hyperbolic, 
t  r  i  1  a  t  e  r  a  t  i.  on  i  or  conventional  navigation  with  chronometer  and  sextant. 
Given  two  lines  of  position  each  with  its  own  standard  deviation,  the 
problem  is  to  determine  the  probability  that  the  meusured  position  is 
within  a  certain  distance  of  the  true  position.  Also  of  interest  is  the 
inverse  problem,  to  determine  the  radius  of  a  circle  around  the  measured 
position  within  which  the  navigator  knows,  with  a  given  probability,  his 
true  position  lies. 

An  explicit  solution  for  the  integral  which  must  be  evaluated  to 
solve  these  two  problems  is  not  possible,  but  the  integrals  can  be 
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evaluated  by  a  digital  computer.  The  mathematical  litoruture  seems  tu 
contain  neither  tubulated  values  of  the  integrals  for  the  nonor thogoru 1 
case  nor  analysis  of  nonortliogonul  bivariate  normal  distributions. 

Material  on  orthogonal  bivuriuto  normal  distributions  is,  however,  com* 
parattvely  ubundunt  (Hefs.  1-8).  In  order  to  avoid  the  necessity  for 
pe r I  or mi ng  many  integrations  by  quadratures,  a  technique  lias  been  found 
which  permits  existing  tubles  of  the  bivuriutu  normal  distribution  to  be 
used  to  solve  the  problems  listed  above.  This  appendix  will  derive  the 
necessary  equations  for  obtaining  cr *  and  cr^ .  The  use  of  these  functions 
as  part  of  Method  2  has  already  been  described.  The  following  appendix 
indicates  how  these  same  functions  may  be  used,  after  calculation  of  some 
auxiliary  functions,  to  obtain  solutions  from  available  tables  other 
than  those  already  illustrated  in  the  body  of  the  memorandum. 

2  .  STATEMENT  Of  THE  PHOBLEM 

Consider  a  general  point  P  whose  position  is  determined  by  measuring 
its  distances,  rj-and  r„  respectively,  from. two  points  of  known  location, 

A/  and  S  (Fig.  B- 1 ) .  M  and  S  might  correspond  to  the  master  and  slave 
stations  of  u  tri lateration  system.  They  might  also  correspond  to  two 
points  ol’  known  location  to  which  the  ranges  r  j  and  r2  are  measured  by  a 
ranging  system  such  as  a  radar,  a  laser,  or  an  optical  rangefinder.  M 
and  S  might  be  points  on  the  surface  of  the  earth  directly  beneuth  two 
stars  that  are  being  used  to  determine  position  by  conventional  navigation 
techniques,  using  sextant  and  chronometer.  The  measurements  of  distance 
are  assumed  to  be  normally  distributed  with  standard  deviations  or1  and  cr z. 

In  fig.  1  concentric  circles  of  ,  r,,  and  rj  -  about  M  and 

r  2  +  ^2'  r2'  anl^  r2  “  a2  ab°ut  5  have  been  drawn  to  suggest  the  normul  dis¬ 
tribution  of  the  two  measurements  of  range.  If  the  sta,  -'  t  '■  deviations  in 
range  are  very  small  in  comparison  with  their  respective  •  •  •■ges,  the  concen¬ 
tric  circles  can  be  assumed  to  be  straight  lines  as  in  fig.  B-2.  In  this 
illustration,  coordinate  axes  Uj  and  are  defined  perpendicular  to  L  j  and 
( 2  respectively,  which  represent  the  true  lines  of  position  corresponds 
to  circles  of  radius  r j  and  r2  in  Fig.  B-l.  Other  lines  of  position  have 
been  added  at  distances  cr  ^  and  cr2  respective  ly  from  1,  and  l2,  With  no 
loss  in  generality,  the  measured  lines  of  position  could  be  assumed  to  be 
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produced  by  conventional  nu v i go t ionu 1  techniques  using  sextant  and  theodolite, 
by  a  tri lateration  system,  or  by  a  hyperbolic  navigation  system.  Random 
errors  in  measurement  may  cause  a  measured  line  of  position  to  bo  displaced 
perpendicular  to  itself,  and  the  probability  that  a  measured  line  of  posi¬ 
tion  l im  will  fall  with  a  zone  of  width  du1  at  a  distance  Uj  from  the  true 
line  of  position  /.  j ,  is  given  by 


P  1^1 

Similarly,  the  probability  that  the  other  measured  line  of  position  i2 
will  fall  within  a  zone  of  width  du2  at  a  distance  u2  from  the  true  line 


of  position 

i2  is  given  by 

p2du2  « 

2 

“2 

1  la  „ 

-  e  du  2 

^2 rrcr2 

(2) 

Notu  that  u, 

and  u j  are  measured 

ulong  the  Uj  and  u2  axus, 

which  in  general 

are  not  orthogonal  except  under  special  conditions,  The  probability  that 
the  measured  position  will  fall  within  an  element  of  area  du  2 ,  du2  centered 
on  the  point  (uj,  u2)  is  given  by 

PJAduldu  j  «  p  lp2du1du2  (3) 

The  probability  that  the  measured  position  falls  within  an  urea  of  any 
size  and  shape,  4,  is  found  by  integrating  the  above  equation  over  that 
area . 


2  a. 


Z^JIC 


(1) 


PA  3  SJ  plp2duldu2  (4) 

Therefore,  in  order  to  find  the  probability  that  the  measured  position 
lies  within  a  distance  r  from  the  true  position,  the  above  integral  must 
be  evaluated  over  a  circular  urea  with  radius  r: 
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1 


du j  du2 


(5) 


2  WjOj 


i(a»r 


The  inverse  problem,  to  find  the  radius,  of  a  circle  centered  on  the  true 
position  within  which  a  meusured  position  has  a  given  probability  of 
fulling,  is  solved  by  fixing  PA  in  the  above  equation  and  solving  for 
r.  Next  we  shull  show  how  existing  tables  can  be  used  to  evaluate  the 
above  nonorthogonul  bivariate  probability  integral, 

3.  EVALUATION  OK  THE  PROBABILITY  INTEGRAL 

In  the  above  integral  Uj  and  u2  are  stochastically  independent 
nonorthogonul  variablesi  The  probability  integral  will  be  solved  by 
first  converting  to  a  new  orthogonal  coordinate  system  with  normally 
distributed  variables.  However,  these  new  orthogonal  variables  will 
not  lie  along  the  major  und  minor  axes  of  the  ellipse.  Hence  the  equa¬ 
tion  of  thu  ellipse,  expressed  in  terms  of  the  new  coordinate  axes  of 
ij  und  will  contain  n  cross-product  term  in  2  ^  z  2 .  Since  such  terms 
are  inconvenient  for  reudy  algebraic  manipulation,  a  second  coordinate 
transformation,  consisting  of  u  simple  rotation,  will  be  employed  to 
remove  this  cross-product  term.  The  final  answers  thus  obtained  will 
be  the  standard  deviations  along  the  mujur  and  minor  axes  of  the 
ellipse, 

Neither  of  the  tuio  coordinate  trails  format  ions  alters  the  shape  of 
the  bivariate  probability  distribution.  After  the  second  transforma¬ 
tion,  the  resulting  probability  integral  can  be  evaluated  by  existing 
tables,  such  us  chose  given  in  Refs.  1  through  \ 

The  first  coordinate  conversion,  illustrated  in  Figs,  B-3,  B-4, 
and  B-5,  is  from  u2  coordinates  to  zL,  j2  coordinates.  The  latter 
two  axes  are  orthogonal  with  z2  coincident  with  Uj.  The  former  two 
axes  are  not  orthogonal,  and  the  and  y2  coordinates  of  the  point 
are  found  by  drawing  two  lines  through  the  point  parallel  to  the  Uj 
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FIG.  B-3  FIRST  COORDINATE  CONVERSION 


und  u2  axes.  Tims  a i  and  u2  coordinates  are  then  found  by  measuring  the 
distance  from  where  these  lines  intersect  the  ul  and  u2  axes  to  the  origin 
of  the  coordinate  system.  By  elementary  analytical  geometry, 


*1  m  21 


kl  +  2j 

u2  -  - 7  (7) 

(1  +  k*f 

where  k  is  a  constant  describing  the  slope  of  line  of  position  l2. 

k  -  t  a  n  Of  ( 8 ) 

where  a  is  the  angle  between  the  two  lines  of  position. 

It  is  necessary  when  converting  from  one  coordinate  system  to  another 
to  replace  the  old  element  of  areu  by  a  new  element  of  area  multiplied  by 
the  Jacobian  of  the  old  variables  witli  respect  to  the  new  variables. 

That  is 


duj  du2 


d  z  j  d  z  2 


66 


By  definition 


D  (  U  J  ,  U  2  ) 

du  j 

du  2 

du  j 

3u 

3  (  Z  1  ,  l  2  ) 

S*2 

(10) 


When  the  indicated  purtiul  derivatives  are  evaluated  hy  partial  differentia 
Lion  of  equations  (6)  and  (7), 


(13) 

The  change  of  variables  has  not  disturbed  the  location  of  any  point,  and 
therefore  the  circular  area  over  which  the  integration  is  to  be  carried 
out  is  unchanged. 

Examination  of  the  exponent  shows  a  cross-product  term  in  z ^  z 2 •  Only 
if  k  were  infinite,  corresponding  to  the  special  case  in  which  the  and 

u2  axes  were  perpendicular,  would  this  term  disappear.  In  order  to  remove 
this  term,  the ' fol lowing  transformation,  which  corresponds  to  a  rotation 
through  an  angle  9 ,  will  be  applied: 

2j  =  x  cos  6  -  y  sin  9  (14) 
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■  T»-8i*?-ia* 

FIG.  B-6  SECOND  AXIS  TRANSFORMATION 


In  the  special  case  when  0"  and  cr2  are  squHl, 

tan  26  *  k  ( 18 ) 

The  slope  of  line  of  position  Z2  with  respect  to  line  of  position  Zj  is 
k ,  and  if  the  angle  included  by  the  two  lines  of  position  is  called  a, 
then 

tan  a  «  k  (19) 

6  =  -  (20) 

2 

Therefore,  when  the  standard  deviations  of  the  distributions  of  the  two 
lines  of  position  are  equal,  the  major  and  minor  axes  of  the  contour  lines 
of  equal  probability  on  the  nonorthogonal  bivariate  probability  distribu¬ 
tion  are  midway  between  the  lines  of  position.  (Compare  Method  1.) 

With  the  cross-product  term  eliminated  by  rotation  through  the  angle 
9  defined  by  equation  (17)  the  desired  probability  is  given  by 


(21) 


From  equation  (16),  the  variances  cr 2  und  cr2  muy  be  found  from 

—  =  - [(ct2  +  o,  +  fc2a2)  cos2  9  +  21(0'?  sin  0  eos  6  +  /;2cr2  sin2  0 ]  (22) 

O'2  0-^(1  *  k2)  1 

und 

2fecr2  sin  0  cos  &  +  fe2o2  cos2  0]  (23) 

Values  of  cr2  and  cr2  are  needed  so  that  equation  (21)  may  be  solved  by  means 
of  existing  tables  of  the  bivariate  normal  distribution  whose  use  will  be 
discussed  in  the  next  appendix.  Once  9  has  been  found  using  equation  (17) 
the  preceding  two  equations  can  be  used  to  calculate  cr2  and  cr2.  It  might 
be  objected  that,  this  is  u  laborious  process.  Fortunately  it  is  possible 
to  obtain  equations  for  cr2  und  cr2  us  simpler  functions  of  cr  j ,  0'2,  und  a, 

By  adding  equations  (22)  and  (23),  thus  obviating  the  need  for  the  auxiliary 
computation  of  the  ungle  0,  we  obtain 


-  ,  - - -  [(cr2  +  cr2  +  fe2(r2)  sin2  8  ~ 

e*  '  crjqgU  +  k2)  ■ 


J_  .  ± 

cr2  cr2 
*  f 

Also  from  equations  (22)  and  (23) 


(24) 


cr2cr2  =  {cTjO2  ( 1  +  k2)}  {[(cr2  T  cr2  +  fe2cr2)  “  4fc2cr{  +  feV4]  sin2  8  cos2  6 
+  [fe3cr4  -  Aa^(a-j  ter2  +  k 2cr2 ) j  sin  9  cos  6  (c.os2  9  “  sin2  9) 

+  k2c72(cr2  +  cr2  +  fc2cr2)(sin4  9  +  cos4  6)}  1  (25) 


From  equation  (17)  it  follows  that 


sin  2  8 


|[21(0'2]2+  [cr2  ( 1  “  k2)  +  cr2(l  +  fe2)]2| 
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1 

2  sin2  a 


[a-2  +  c r\ 


(o2  +  ir\)  -  4(«in2  a)ujo-2l 


With  the  aid  of  tables  of  squares  anti  squares  of  sines  of  unglos  (He  I',  10, 

for  example),  these  equations  muy  be  solved  to  determine  numericul  values, 

It  is  to  be  noted  that  the  equations  are  in  terms  of  variances;  to  obtuin 
stundurd  deviations  as  desired  for  some  calculations  square  roots  must  be 
taken.  It  is  also  to  be  rioted  thut  the  formulas  for  the  two  desired  func¬ 
tions  differ  only  in  the  sign  before  the  radical  sign.  Thus  in  numerical 
calculation  the  solution  of  one  expression  provides  all  the  numbers  required 
for  the  second. 

For  the  special  case  when  cr1  =  a 2  equations  (36)  and  (37)  may  be  greatly 
simplified  as  follows; 


2  sin2  a 


(2cr2  +  / 4cr*  -  4  s  in  2  tttr  * } 


{; 2<-r 2  +  2cr2  /l  -  sin2  d} 


2  sin2  a 


2  sin2  a 


(2 cr1  tl  +  cos  aj  ) 


Since 


1  +  cos  a 


sin  a 


1  +  cos  a 


1  +  cos  ci 


sin  a  »  sin  2(a/2) 


2  sin  — 

2 


i  v 


and  s  imi  lui'  ly 


o/2 


(46) 


The  ratio  &y/vx  (s  c  of  Method  2,  Section  2)  is  useful  in  further  conver¬ 
sion  of  the  error  ellipse  to  circles  of  equivalent  probability  and  may  be 
readily  seen  to  be 


tun  1/2  a 


when  cr1  =  ct2 


(47) 


4  .  ALTERNATIVE  FORMULAS 


Reference  9  lists  the  following  formulas  for  cr*  and  cr  ^ 

1 


2 


(c2  +  o-j )  esc2  a  +  0 1  rot  a  esc  20 


r 1  =  —  (a2  +  cr2  )  CSC2  a  -  cr  cot  a  CSC  20 

y  2  1  2  1 


tan  20  - 


cr2  sin  2a 


cr2  COS  2a  +  or2 


6  is  the  angle  between  cr1  and  <J x  us  in  the  foregoing  analysis, 

Although  the  expressions  of  equations  (36)  and  (37)  initially  appear 
lengthy,  experience  Iihs  shown  thut  they  are  well  suited  to  computation. 
They  have  the  udvnutuge,  too,  not  possessed  by  the  alternative  formulas, 
of  not  requiting  uuxiliary  computation  of  additional  functions  before  the 
given  values  may  be  entered  into  the  formulas.  Equations  (36)  and  (37) 
utilize  only  the  basic  input  data— the  errors  along  the  given  axes  and 
the  angle  between  these  axes.  Avuilable  handbooks  of  mathematical  tables 
(Ref.  10,  for  example)  contain  tables  of  squares  and  also  of  the  squares 
of  sines  of  angles  so  thut  the  numerical  values  needed  may  be  readily 
obtained.  Although  a  desk  calculator  is  handy,  slide  rule  computations 
are  not  difficult  and  are  of  sufficient  accuracy  for  systems  evaluation. 
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BIVARIATE  PROBABILITY  DISTRIBUTION 
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APPENDIX  C 


EXAMPLES  OF  THE  USE  OF  TABLES  OF  THE 
BIVARIATE  PROBABILITY  DISTRIBUTION 


1.  INTRODUCTION 

Reference  1  contains  the  most  extensive  set  of  tables  applicable  to 
bivariate  normal  distributions.  The  use  of  these  tables  requires  the 
calculation  of  certain  auxiliary  functions,  which  are  described  in  this 
section.  A  number  of  problems  illustrating  the  use  of  these  functions 
and  tables  follow, 

Let  Pr(ajUj  +  OjUj  -  t}  denote  the  probability  that  a  randomly  chosen 
point  will  fall  inside  an  ellipse  whose  equation  is 

+  a  2u\  »  t  (1) 

where  and  a2  are  stochastically  independent  normally  distributed  variables 
with  means  equal  to  zero  and  standard  deviations  equal  to  one.  By  defi¬ 
nition, 


a,  +  a, 


«  1 


then 

Pjfoj.aj.'t)  =  P r(a jU j  +  a 2u j  =:  t} 


and 


P  2  ^  B  1 '  2  '  ^  ^ 


1 


2W„  lCTu  2 


a 


2  \ 
u  ]  ■* 


1  2  u  2 


(fu^O iuj 


(2) 

(3) 


(4) 


where  cr  »  cr  =>  1, 

u  i  u  2 

Tables  of  PJ(.al,ai ;  t)  are  given  in  Refs,  1,  2  and  3.  Reference  1 
gives  120  values  of  P2  (a^  a2;  t )  to  four  significant,  figures,  with 
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Preceding  page  blank 


t  15  0, 1(0.1)  1.0(0, 5)  2. 0(1.0)  5.0.  (This  inruns  that,  t  has  intervals  of 
0.1  be  tween  0.1  und  1.0,  and  0.5  between  l.Ound  2,0  etc.)  and  with 
s2  a  0,5,  0.6,  0,7,  0.8,  0.9,  0.95,  0.99  und  1.00.  Reference  2  gives 
149  values  of  P2 (a1# a2;  t  )  to  five  significant  figures,  with  t  ~  0.1(0. 1) 
1.0(0, 5)  2. 0(1.0)  5.0  und  with  a2  -  0.5,  0.6,  2/3,  0.7,  0.8,  0.875,  0.9, 
0.95  and  0,99.  Another  table  in  Ref.  3  gives  150  valuos  of  the  inverse 
of  u  quadratic  form  in  two  dimensions,  i.e.,  the  value  of  t  which  causes 
P^(alia2'  t)  to  huve  u  certain  value.  The  values  of  u2  are  the  same  as 
in  the  preceding  table,  and  P  =  0.05(0.05)  0.30(0.10)  0.70(0.05)  0.9b. 
Reference  3  gives  2881  valuos  of  (u n2; t )  to  eight  significant  figures 
with  t  =  0.005  (0.005  )  0.10(0.01  )  1. 00(0. 02  )■  2. 5  0(0. 025  )  3.50(0.05) 
5.00(0.25)  6.00(0,50)  7.00(1.00)  10.00  and  a?  =  0.5,  0,6,  2/3,  0.7,  0.75, 
0.8,  0.875,  0.9,  0.95,  0.99  and  1.00. 

Thenc  tables  can  still  be  used  even  if  the  two  variables  do  not  have 
unit  standard  deviations.  Consider  two  normally  distributed  variables 
and  *2  with  standard  deviations  and  cr •  respectively.  If  the  bivariate 
probability  integral  is  evaluated  over  t lie  area  inside  an  ellipse 


6 1*1  + 


.  x2 
>2*2 


(5) 


then 


Pr^bi> 


+  6  2^ |  =  r2} 


1 


2  rrer 


b 


*l+fc2x2J‘r 


(6) 


In  order  that 


^  i  x  i  +  b2xl  = 


M^u*  +  a2u2  =  t) 


it  is  necessary  that 


ai«l 


b2X  2 


a2Ul 


t 


p  2  ( a  i  < a  2 ;  t )  ( 7 ) 


(8) 


where  .  a  constant. 
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Since  x  j ,  .r2,  ii|,  unii  u2  are  nil  normally  distributed 


and 


1 1 

z 


u  1 


u: 


au  °-«2 


Since  a  .  and  cr  equal  one, 

u  1  .  u  2  * 


and 


u l 


*  2 


When  these  values  are  substituted  into  equation  (10) 


blul\ 


b2ali 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


and 


(15) 


Since  by  definition  the  sum  of  a1  and  u2  is  one, 

biatl  b 2°  j 2 


(16) 


or 


6l°1l  * 


(17) 
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There  Lore 


P2 (dj, a2; t ) 


27 W',  1°',  2 


*  1  *  2/ (/*  [f/jf  2 


(18) 


b  1*  l+b2,2[ 


where  Oj,  a2  and  t  us  functions  of  b , ,  b2,  O'2  ,  ,  ^2und  r2  are  given  by 
equations  (13),  (U),  (15),  and  (16).  These  equations  disagree  with  the 

ones  given  in  Ref.  2;  as  the  above  derivation  shows,  Ref.  2  is  wrong. 

When  equation  (20)  is  substituted  into  equation  (21)  of  Appendix  B, 


K  1°-,  2 


o2a-  ( 1  +  k2 ) 

1  i! 


P2 («! , «2;  t  ) 


(19) 


But  from  equation  (31),  Appendix  B, 


(1  +  feM 


}A 


c  r  ,CT  ,, 

t  1  X  2 


^2 


(20) 


so 


^  2  (  ^  i  <  ^  ^ 


(21) 


Therefore,  in  the  general  case,  PA  ,  the  probability  of  being  within  an 
ellipse,  b^2  +  b2*2  ■  r2,  can  be  read  directly  from  tables  of  P^a^a^.t) 
after  ax,  o2,  and  t  have  been  calculated  means  of  equations  (36)  and  (37) 
of  Appendix  B  and  (13),  (14),  (15),  and  (16).  In  the  special  case,  with 
which  this  memorandum  is  concerned,  finding  the  probability  that  the  true 
position  is  within  a  distance  r  of  the  indicated  position,  the  calculations 
are  still  simpler.  The  ellipse  becomes  a  circle,  bj*2  +  b2*22  =  r2,  and 


By 


b,  -  b2  =  1 

substitution  using  equations  (35)  of  Appendix  B  and 


(22) 


+  crt 


sin  a 


(23) 
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r.qua  Lions  (13)  to  (15)  simplify  to 


4  ( s  i  a  2 

a)  a 

2  2‘ 

CT 

1  2 

<«■? 

9 

T  (J  * 

)2  _ 

4  ( s  i  a  2 

a  )o 

\ul 

K 

)2 

r  2  sin2  a 

a2  +  cr2 

2  .  PROBLEMS 

Table  C-l  summarizes  the  values  of  the.  given  quantities  in  the  dif¬ 
ferent  problems  und  the  resulting  probubi 1 . t.os  P2(a1,a2;t).  Where  possible 
the  problems  have  been  cheeked  by  use  of  different  tables.  As  a  convenience 
to  the  reader,  Tables  C-2,  C-3,  and  C-4,  abstracted  from  Ref.  2,  are 
included  at  the  end  of  this  appendix  su  that  the  results  presented  here 
cun  be  verified.  Table  C-2  is  n  brief  table  of  P2(a1,a2;(). 


Table  C.-1 

SUMMARY  OF  PR0B1.KMS 


PHOBI.EMS 

a 

CT2 

*2 

r 

^2<°1*  °2’  * 1 

1 

90° 

1 

1 

1 

1 

i 

.0.39347 

2 

90° 

1 

1 

1 

1 

2 

0.86466 

3 

90° 

2 

2 

1 

1 

2 

0.39347 

4 

90° 

1 

1 

1 

2 

2 

0.74244 

5 

90° 

1 

2 

1 

1 

1 

0.21529 

6 

90° 

1 

2 

1 

1 

2 

0.59009 

7 

90° 

1 

2 

4 

1 

2 

0.39347 

8 

90“ 

1 

2 

1 

4 

2 

0.32623 

9 

60“ 

1 

1 

1 

1 

1 

0.34230 

10 

60“ 

1 

1 

1 

1 

/273 

0.24601 

11 

60“ 

1 

1 

1 

1 

VI 

0.55620 

12 

60° 

1 

1 

1/2 

3/2 

1 

0.39347 

13 

60° 

l 

1 

3/2 

1/2 

1 

0.34945 

14 

65°  700  ft 

900  ft 

1 

1 

1520  ft 

0.77624 

15 

65°  700  ft 

900  ft 

1 

t 

3040  ft 

0.99404 

The  more  exLensive  Luiilo  of  P2  (n  t ,  a2 ;  t  )  in  Hef,  3  wus  actually  used  in 
working  these  problems  because  it  hus  more  entries,  und  the  interval  for 
interpolation  is  smaller.  Table  2  can  be  used  to  integrate  (1)  a  circular 
probability  distrioution  over  a  circular  urea,  (2)  a  circular  probability 
distribution  over  an  elliptical  area,  (3)  an  elliptical  probability  distribu 
t.ion  over  u  circular  area,  or  (4)  an  elliptical  prooubility  distribution 
over  an  elliptical  urea.  A  bivariate  normal  distribution  such  us 


is  culled  a  circular  distribution  when  -  cr 2  because  the  contours  of 

equal  probability  density  are  circular.  If  cr  l  f  the  lines  of  equal 

probability  density. ure  elliptical  und  the  bivariate  distribution  is 
called  an  elliptical  distribution.  Table  C.-2  can  be  used  to  integrate 
a  circular  probability  over  a  circular  area,  and  Table  C- 3  can  be  used 
(1)  to  integrate  a  circular  probability  distribution  over  a  circular  area, 
or  (2)  at  the  expense  of  entering  the  table  twice,  to  integrate  a  circular 
probability  distribution  over  an  elliptical  area.  Since  the  latter  table 
occupies  several  pages,  on  L y  those  pages  needed  to  work  the  problems  in 
this  memorandum  have  been  included.  As  a  check,  several  problems  have 
been  solved  using  Tables  C-3  or  C.-4  or  both  as  a  check  an  the  answer 
obtained  from  Table  C-2.  Unfortunately  not  ull  the  problems  could  be  so 
checked,  because  of  the  inherent  limitations  of  Tables  C-3  and  C-4. 

Prob  lem  ( 1) 

Given 

a  =  90°,  crt  =  ut]  =  1,  a 2  =  <tx2  =  1,  6j  *  1,  b2  -  1,  r  -  1 

Find 

PjUj,  a2;  t ) 

Since  the  lines  of  position  are  perpendicular  [a  =  90°)  and  their  standard 
deviations  are  equal,  the  probability  distribution  is  circular.  The  area 
over  which  the  integral  is  taken  is  likewise  circular,  (6j*^  +  b2*2  =  rZ 
or  x j  +  *2  =1)  and  the  area  has  a  radius  equal  to  one  standard  deviation. 


82 


Table  C-2 

DISTRIBUTION  OF  A  QUADRATIC  FORM  IN  TWO  DIMKNSIGNS 


*2*  'i 

0.1 

0.2 

0.3 

0.4 

0.3 

,3  ,  .3 

0. 093U 

0.10127 

0.23910 

0.32900 

0,39347 

.3  ,  .4 

0.09093 

0.10429 

0.20304 

0.33403 

0.39009 

2/3,  1/3 

0.10033 

0.19003 

0.27033 

0.34222 

0.40004 

,7  ,  .3 

0.10200 

0.19432 

0.27300 

0.34013 

0.41270 

>73  ,  .23 

0.10014 

0,20299 

0.20037 

0.33902 

0.42440 

>1  ,  .2 

0.11301 

0.21329 

0.30112 

0.37330 

0.44023 

.373,  .123 

0.13340 

0.24401 

0.33434 

0.40040 

0.47117 

,9  .  .1 

0,14000 

0.23941 

0.34943 

0.42237 

0.40313 

,93  ,  .03 

0.10130 

0.30010 

0.30301 

0.43200 

0.30594 

,99  ,  .01 

0.23300 

0.31030 

0.41130 

0.40900 

0.31020 

0-0 

0.7 

Q.» 

-AiA- 

,3  ,  .3 

0.43119 

0.30341 

0.33007 

0.39343 

0.43212 

,1  ,  .4 

0.43305 

0.30797 

0.33300 

0.39740 

0.03379 

2/3,  1/3 

0.44441 

0.31023 

0.30179 

0.00402 

0.04223 

>7  .  .3 

0.47040 

0.52203 

0.30019 

0.00932 

0.04430 

,73  ,  .25 

0.40200 

0.33294 

0.37014 

0.01039 

0.43429 

,9  .  .2 

0.49070 

0.34040 

0.39009 

0.02072 

0.00297 

,075,  .123 

0.32433 

0.57011 

0.00905 

0.04407 

0.07340 

,9  ,  .1 

0.33423 

0.37793 

0.01307 

0.04910 

0.07040 

,93  ,  .03 

0.33133 

0.39040 

0.02400 

0,03401 

0.00192 

,99  ,  .01 

0.33902 

0.39013 

0,42027 

0.03092 

0.40200 

■  -1.8- 

.  Ml.. 

■■  \A... 

-JUL. 

3  ,  .3 

0.77007 

0.00400 

0.93021 

0.90140 

0.99320 

4  ,  .4 

0.77049 

.  0.00401 

0.94071 

0.90023 

0.99227 

2/3,  1/3 

0.70100 

0.00424 

0.94000 

0.97773 

0.99037 

7  ,  .3 

0.70202 

0.00379 

0.*4412 

0.97400 

0.90940 

73  ,  .23 

0.70491 

0.00233 

0.94C04 

0.97310 

0.90732 

,0  .  .2 

0.70004 

0.00030 

0.93033 

0.90904 

0.90331 

073,  .123 

0.70070 

0.03300 

0.92944 

0.90430 

0.90100 

,9  ,  .1 

0.70301 

0.03274 

0.92093 

0.90244 

0.90020 

,93  ,  .03 

0.70290 

0.04704 

0.92107 

0.93131 

0.97733 

99  ,  .01 

0.70009 

0.04374 

0.91770 

0.93331 
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Tttbltj  C-3 

CRITICAL  VALIJKS  KOR  THK  CIRCULAR  NORMAL  DISTRIBUTION 


JL_ 

I 

P 

» 

J> 

a 

O.Ol 

0.1414 

0.44 

1.1101 

0.91 

2.1943 

0.02 

0.2010 

0.47 

1.1248 

0.92 

2.2473 

0.03 

0.2448 

0.48 

1.1U4 

0.9) 

2.3042 

0.04 

0.2437 

0.49 

1.1403 

0.94 

2,3721 

0.0) 

0.3203 

0.30 

1.1774 

0.93 

2.4477 

0.04 

0.3314 

0.31 

1.1944 

0.94 

2.3373 

0.07 

0.3410 

0.32 

1.2114 

0.97 

2.6U2 

o.oa 

0.4044 

0.33 

1.2244 

0.94 

2.7971 

0.09 

0.4343 

0.34 

1.2442 

0.99 

3.0349 

0.10 

0.4390 

0.33 

1.2637 

0.991 

3.0494 

0.11 

0.MZ4 

0.34 

1.2814 

0.992 

3.1073 

0.12 

0.3034 

0.37 

1.2992 

0.993 

3.1302 

0.13 

0.3274 

0.34 

1.3172 

0.994 

3.1947 

0.14 

0.3492 

0.39 

1.3334 

C.99S 

3.13)2 

0.13 

0.3701 

0.40 

1.35)7 

0.994 

3.3231 

0.14 

0.3903 

0.41 

1.3723 

0.997 

3.4084 

0.17 

0.4103 

0.42 

1.3911 

0.994 

3.32)3 

0.14 

0.4)00 

0.4) 

1.4101 

0.999 

3.7149 

O.lt 

0.4492 

0.44 

1.U94 

0.9991 

3.7432 

0.20 

0.4440 

0.43 

1.4490 

0.9992 

3.7743 

0.21 

0.4444 

0.44 

1.4449 

0.9993 

3.8117 

0.22 

0.7049 

0.47 

1.4891 

0.9994 

3.8319 

0.23 

0.7230 

0.44 

1.3094 

0.9993 

3.8989 

0.24 

0.7409 

0.49 

1.3303 

0.9994 

3.9338 

0.23 

0.7343 

0.70 

1.3314 

0.9997 

4.0278 

0.24 

0.7740 

0.71 

1.3733 

0.9998 

4.1273 

0.27 

0.79)4 

0.72 

1.3934 

0.9999 

4.2919 

0.24 

0.4104 

0.73 

1.4182 

0.99993 

4.4303 

0.29 

0.4274 

0.74 

1.4414 

0.99999 

4.7983 

0.30 

0.4444 

0.73 

1.4431 

0.999993 

4.9409 

0.31 

0.4413 

0.74 

1.4494 

0.999999 

3.2343 

0.32 

0,4713 

0.77 

1.7143 

0.9999993 

5.3444 

0.33 

0.4930 

0.74 

1.7402 

0.9999999 

3.4777 

0.34 

0.9114 

0.79 

1,7447 

0.99999993 

5.7983 

0.33 

0.9242 

0.40 

1.7941 

0.99999999 

4.0497 

0.34 

0.9444 

0,41 

1.8223 

0.999999993 

4.1429 

0.37 

0.9413 

0.42 

1.8319 

0.999999999 

4.U79 

0.34 

0.9774 

O.U 

1.4423 

0.39 

Q.99U 

0.44 

1.9143 

0.40 

1.0104 

0.43 

1.9479 

0,41 

1.0273 

0.84 

1.9830 

O.U 

1.0U4 

0.47 

2.0200 

O.U 

1.0403 

O.U 

2.039) 

0.44 

1.0749 

0.49 

2.1011 

0.4) 

1.0933 

0.90 

2.1440 

NOT  REPRODUCIBLE 
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NOT  REPRODUCIBLE 

I'll  I,  I  r  C-4 

Ol'KSKT  CUV.IJi  1'HOHAHII.ITIIiS  I'OH  Till',  C.IHClil.AH  MOMMA).  iilSTHIIlDTION 


(r«  -  Dl/ff 


-  V?, 

+  - 

(l)(l)  +  (1)0)  ■  2 

1 a  «  1 

(1)(1) 

1 

c 

2 

2 

(1)(1) 

1 

From  Table  C-2,  P2(a,,a2;t)  ■  0.39347.  ThiB  problem  can  be  checked  by 
Table  C-.3.  This  table  gives  corresponding  values  of  P  and  fi  where  P  is 
tlie  probability  under  the  circular  normal  distribution  for  a  given  value 
of  B,  und  B  =  r/a  where  r  is  the  radius  over  which  the  circular  normal 
probability  distribution  is  integrated,  and  cr  is  the  standard  deviation 
of  the  eirculur  normal  distribution  measured  along  any  axis.  For  the 
values  given  above,  B  =  r'/cr  -  1/1  =  1.000,  From  Table  C- 3 , 


P  B 


0.39 

0.9943 

0.40 

1  010B 

0.01 

0.016D 

By  interpolation, 

/0. Q057\ 

pn-i  ooo  -  0.39  +  -  (0.01)  »  0.39346 

B -1'000  \0. 01 65/ 

This  checks  the  value  obtained  above. 

A  second  check  is  possible  using  Table  C-4.  This  table  gives  q(rrf/cr,  D/cr) 
which  is  the  probability  that  a  randomly  chosen  point  will  fall  outside  a 
circle  of  radius  r d  with  center  at  (£J , 0 )  when  the  underlying  probability 
distribution  is  a  circular  normal  distribution  centered  at  the  origin  and 
with  its  two  standard  deviations  equal  to  cr.  Then 


(27) 
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will'  IT 


llio  ii  rou  ovi'i'  which  t-his  integral  is  tukcti, 


(x  -  I))'2  +  y2  =  r'1 


(28) 


(29) 


is  tin'  equation  of  it  circle  of  radius  r  d  centered  at  (D,  0).  If  0  is  set 
equal  to  zero,  r j  -  1  anil  the  standard  deviation  of  the  circular  normal 
distribution  set  equal  to  1,  then  q{rd/o\,  I)  ]o )  can  be  read  from  Table  C-4 
and  the  desired  probability  found  for  equation  (27).  The  table  of 
ijti1  ,/u,  D/a)  uses  (r^-M/cr  and  O/o  ns  arguments. 


(1)  ~  (0) 
1 


*  1.000 


D 
'  a 


0 

-  =  0.000 
l 


From  Table  C-3,  i/(l,0)  =  0.607,  so  p  ( 1 , 0 )  =  1-  9(1,0)  =  1  -  0.607  =  0.393 
which  checks  the  two  vulues  previously  found. 


Problem  (2) 
Qi  veu 


a  =  90°,  ox  =  crj  l  =  1,  cr2  =  o  j2  =  1,  6, 


1,  b2  -  1,  r  -  2 


F  i  ii  d 


P  2  (  11  j  »  Q  2 ; f  ) 


This  problem  uses  the  sume  circular  probability  distribution  as  was  used 
in  problem  (l),  but  the  distribution  is  integrated  over  n  circle  whose 
radius  is  twice  the  radius  of  the  circle  used  in  problem  (1).  Consequently 
the  probability  found  in  this  problem  should  be  greater  thun  the  probabil¬ 
ity  found  in  problem  (1). 


C 


a 


l 


6,^,  T  bitrU  "  (1)(1)  +  (1X1)  ■  2 


_  (1)(1)  _  l 

C  2  2 


07 


b  ,,u  2 

1  X  1 


(1)(1) 

2 

=  2 


1 

2 


From  Table  ('-2,  P^  (a  (  ,  «2  ;  t  )  =>  P^  (0. 5 , 0.  5;  2  )  B  0,86466 

Hi  i  s  answer  cun  Ijo  oliookoil  by  Table  C-3.  H  -  r/o'  =  2/1  ~  2.  From  Table  C-3, 


P  B 

0.86  1.0830 

0,87  2.0200 

0.01  0.0370 


By  in  terpol  u  1, i on  , 


B-2 


0.86  + 


0.  017 
0.037 


(0.01) 


0.8646 


A  second  check  is  possible  using  Table  C-4.  From  the  given  quantities, 


rd  -  D  (2)  -  (0) 

a  =  1 


2.000 


D 


0.000 


From  the  table 


and 


q(2.000,  0.000)  =  0.135 


1.000  0.135 


0.865 


Problem  ( .? ) 
Given 


a  =  90°, 


2,  cr. 


x  2 


2,  6,  -  1,  6, 


1, 


=  2. 
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I'i  ml 

l  2  (  H  j  #  2  1  ^  ) 

This  problem  is  similar  to  the  first  problem  in  tliut  tin:  probability  dis¬ 
tribution  is  circular,  a l though  in  this  problem  the  circular  probability 
il  i  s  tr  i  but  i  on  is  wider  and  its  maximum  is  lower  than  tin*  circular  distribu¬ 


tion  in  problem  (l).  line  a use  the  area  over  which  the  integral  in  tukon 
is  increased  in  the  same  proportion  us  the.  standard  deviation,  the  prob- 


u  b  i  1 

1  i t )  in  this  p rob i cm  i > 

;  the  same  as 

the  probnb i  1  i ty 

found  in  problem 

C  -  b ,a* , 

( 1 ) (4  )  +  (1)(4) 

=  8 

(1)(4) 

1 

1,1  5  c 

8 

2 

( 1 ) (4  ) 

1 

a  -c 

8 

2 

r  2 

4  _  l 

l  -  3 

a 

8  ’  2 

Front 

Table  C-2,  P2(aj,a,;t)  =  P2(0.5,0 

.5:0.5)  -  0.39347. 

This  answer 

can 

be  checked  by  Table  C- 

3.  R  *  r/cr  ■ 

=  2/2  =  1.  For  B 

*  1, 

/’  = 

0.39346  as  shown  ubove 

lor  problem 

( 1 ) .  A  check  by 

means  of 

Tuble  C-4  is  also  possible. 
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I’rtth  /cm  f  4  ) 

Give  n 

nt  »  90°,  a-j  =  ax  (  »  l,  Oj  =  ct2  «  1,  lij  =  1,  bj  s  2,  r  =  2 

Kind 

P2  (° i <  u2 ; 1 ) 


l'li  i  s  [i  roll  1  cm  uses  the  same  circular  probability  distribution  as  was  used 
in  problem  (1),  but  the  urea  is  now  evaluated  over  an  ellipse  whose  semi- 
mujor  axis,  2  units  long,  lies  along  the  Xj  axis  and  whose  semiminor  axis 
1  unit  long  lies  along  the  x2  axis.  Since  the  same  circular  probability 
distribution  is  used  in  problems, (1),  (2)  and  (4)  and  the  area  of  Integra 
lion  in  problem  (4)  is  intermediate  between  the  areas  of  integration  used 
in  problems  (1)  and  (2),  the  probability  in  this  problem  should  fall  be¬ 
tween  the  probabilities  found  in  problems  (1)  and  (2). 


G  *  bj  0-*,+  H2<j\2  =  ( 1 )  ( 1 )  +  (2 )  (1 ) 

6ia»l  _  ( 1 H  l )  _  _1 

' 1  C  "  3  3 

b2atx2  _  (  2  )  (  1  )  _  2 


/I  2  4  \ 

-  P2h.  ~J  =  0.74244 


3 


Table  C-3  cannot  be  used  because  it  is  applicable  only  when  the  area  of 
integration  is  circular.  Toole  C- 4  can  be  used  to  handle  integration  of 
u  circular  normal  distribution  over  an  elliptical  ureu.  Let 
Prl6j*j  +  b  2x22  =  r2}  denote  the  probability  that  a  randomly  chosen  point 
will  lie  within  an  elliptical  area  bjX2  +  6 2 x |  =  r2  when  *j  and  x2  are 
normally  distributed  with  unitary  standard  deviations,  which  is  the  same 
us  the  integral  of  a  circular  normal  distribution  .with  the  same  standard 
deviation  integrated  over  the  area  defined  by  that  ellipse.  Then 


Pr{6jX2  +  bjX2  =  r2}  =  q(B,  4)  -  g(4,  B) 


90 


win1  ri' 


2' 

'—I  + 

it., 


1 


\\  <4 

+  I  - 
2 


1.707 


anil 


,2\ 


l; 


4\^ 
2 


0.293 


(f(Bi  A)  -  q 


■ 


In  order  to  find  q(B,  A)  Table  C-4  must  be  entered  usi-g  the  values 
(rj  ~  D/cr  and  D/cr  for  the  arguments. 


q{B,  4)  =  q 


fe.7) 


r.n  “  D\ 


0.293  -  1.707  =  -1.414 


By  interpolation, 


—  =  1.707 

a 


q(B,  A)  =  0.990 


r  d  2  ^2 

q  (A ,  B)  =  q  I - ,  —I  =  g(1.707,  0.293) 

a-  a  I 


'd2  ~  ^2 


1.707  “  0.293  =  1.414 


D2 

—  =  0.293 

cr 


By  interpolation, 


q(A,  B)  -  0.248 
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I'i  null  v 


Wi  +  *.* 


2  < 
2-'  2  = 


</(B,  4)  -  ,;(4,  B) 

0.990  -  0.248  =  0.  742 


Tli  is  chucks  the  previously  obtuiucd  vuluu  of  /’2(a1<a2;t) 


/’/  oh  l  fm  ( ) 
Ui  Vc  I) 


a  *  90 ”, 

CTi  "  = 

>.  ^2  “  cr,  j  =2,  l 

>1  =  1,  62  =  1,  r  «  1 

F  i  n  d 

P2(W 

t ) 

lloouiisu  the  standard  deviations  of  the  bivariate  normal  probability  dis¬ 
tribution  are  unequal,  the  contour  lines  of  the  distribution  are  ellipses 
with  their  mujor  axes  lying  along  the  *2  BXls>  The  integral  of  the  'ellip¬ 
tical  probability  distribution  will  be  found  over  u  circle  whose  radius  is 
one . 


bi°,2i  ■  bi(TU  *  umu  A  (1)(4)  B  5 
( 1) ( 1 ) 


h 

"  1  x  l 


b  crJ 

2 x  2 


5 

(1X4) 


0.2 


0.8 


0.2 


/’jUj,  a2;  t )  =  P2  (0. 2 , 0. 8;  0. 2  }  -  0.21529 


P rob l urn  (6 ) 

Given 

a  -  90°,  cr1  =  o-x  J  =  1,  cr2  »  crj2  =  2,  b  j  =  1 ,  6  s  »  1,  r  »  2 


92 


V  i  n  d 


P2  (<ij.  a2\  t) 

This  is  the  same  as  the  preceding  problem,  except  that  the  radius  of  the 
circle  over  which  the  elliptical  distribution  is  to  be  integrated  has  been 
doubled. 


C  =  +  ba<r\2  -  ( 1 )  ( 1 )  +  (1 )  (4  )  «  5 


b.cr 


1“*  1 


0.2 


a 


2 


b2a 


2 

<2 


c 


0.8 


P2(alta2-,t)  -  P2(  0.2,0.810.8) 


0.59009 


Problem  (7) 

Given 

a  =  90°,  CTj  =  <Tti  =  1,  cr2  =  C7x2  =  2,  b1  «  4,  b2  =  l'  r  =2 

Find 

P2(a j, a2 ; t  ) 


The  elliptical  distribution  is  the  same  as  in  the  preceding  two  problems, 
but  now  the  integral  of  the  elliptical  distribution  is  to  be  integrated 
over  an  ellipse  whose  major  axis  is  parallel  to  the  major  axes  of  the 
ellipses  of  equal  probability  density  on  the  elliptical  distribution. 


C  -  6^ 

+  62o-22  = 

(4X1)  +  ( 1 )  (4) 

(4 )  (1 ) 

-  o.s 

A  1  " 

1  C 

8 

b2al2 

(1X4) 

-  0.5 

c 

8 
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i 


c 

/J2(a1(a2;t)  -  />2  ( 0 . 5 , 0 . 5  i  0. 5  )  =  0.39347 

Problem  (S) 

Given 

a  =  90°,  3-j  =  (7i2  =  1,  cfj  =  (7iS  =  2,  6j  =  1,  62  =  4,  r  =  2 

Find 

P  2  (  n  1  *  a  2  '  i  ^ 

This  problem  uses  the  sane  e 1 1  ipt ica 1  probe bi 1 i ty  distribution  and  the 
sume  elliptical  area  over  which  the  elliptical  distribution  is  to  be 
integrated,  except  that  in  this  problem  the  elliptical  area  has  been 
turned  90°  so  that  its  major  axis  is  parallel  to  the  minor  axes  of  the 
coutburs  of  equal  probability  density  of  the  elliptical  probability 
distribution. 


c  = 

+  b2cr2t2  - 

(1)U)  +  (4 )  (4 )  ■  17 

a  1  3 

( 1 )  ( 1  ) 

=  0.05B82 

c 

17 

a  2 

(4) (4) 

16 

=  —  =  0.9412 

17 

c 

17 

t 

it 

4 

—  =  0.2353 

17 

P2(alta2\t)  =  P2(0. 05882, 0.9412;0.2353)  -  0.32623 

Problem  (9) 

Given 

d  -  60°,  =  1,  cr2  =  1,  bj  =  1,  62  =  1,  r  "  1 

Find 

Pi(al,a2;t) 


0.5 
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Ill  nil  of  the  previous  problems  the  lines  of  position  have  been  assumed 
to  intersect  ut  right  angles.  In  this  and  the  remaining  problems  in  this 
memorandum,  angles  of  intersection  other  than  90°  will  be  used.  The  area 
over  which  the  elliptical  probability  distribution  is  to  be  integrated 
is  a  circle  of  unitary  radius  centered  on  the  origin,  so  equations  (24) 
to  (26)  can  be  used. 


5  <l* 


4  (sin2  <X)a-*cr22 


(< y\  +  cr2)' 


JO 


1  + 


1  - 


1  - 


1  xl 


4(4)  (1X1) 


(1  +  l)z  J 


1  + 


4(sin2  0f)cr2cr 


,  [A 


t  = 


2  ■  2 
r  sin  a 


erf  +  cri 


(cr2  +  cr2) 


[;]'}•  i 

1  -  eft  •  i 


1  +  1 


8 


P2(a,.a2;t)  *  Ps  { 0 . 75 , 0 . 25 ;  0 . 3  75 )  -  0.34230 


It  is  also  possible  to  solve  this  problem  by  first  solving  for  c r*  j  and 
ct22  and  then  using  equations  (14)  to  (17),  as  in  the  previous  problems. 


a 


2 

X  1 


(7 


2 

X  2 


{ O' \  +  cr2  +  [(cr2  +  cr2  )  - 


4  (sin2  oOo-Jcr2]^} 


- - - {cr2  +  cr2  “  [(cr2  t  cr2  ) 2  -  4(sin2  0t  )cr2cr^  ]  '*}■ 

2  sin2  a 


2 

3 


From  these  values,  the  contours  of  equal  probability  density  are  ellipses 
satisfying  the  equation  (*2/2)  +  ( 3 * | / 2 )  =  k2.  The  ellipse  with  fen=  I 
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Problem  (10) 

Given 

a  =  60°,  cr1  =  1,  a2  »  1,  b:  »  1,  6S 

Find 


'I'hc  elliptical  probability  distribution  in  this  problem  is  the  same  ns 
in  the  preceding  problem,  and  the  only  difference  is  that  the  radius  of 
the  circle  over  which  the  elliptical  probability  distribution  is  to  be 
integrated  is  reduced  to  0,8165.  This  circle  is  also  shown  in  Fig.  C-l. 
Therefore  a1  and  a2  are  the  same  as  in  the  preceding  problem. 


t 


a j  =  0.75 

a  2  =  0..25 

r>  sin1  «  (D(?)  „  .. 

-  ,  -  .  0.25 

a2  +  o-2  1  +  1 

■  cr j 


P8(«l#a,st)  ■  P,(0. 75, 0.25; 0.25)  -  0.24601 


Prob  lem  (11) 

Gi  ven 

a  =  60°,  crl  ■  1,  o-j  ■  1,  6j  ■  1,  b2  *  1,  r  *  ^2 

Find 

1*2  (  ai*  a2  ' 


This  problem  is  the  same  as  the  preceding  two,  except  that  the  radius  of 
the  area  over  which  the  elliptical  probability  distribution  is  to  be 
integrated  is  now  1.414.  This  circle  is  also  shown  in  Fig.  C-l. 


°l  ■  0.75 

a2  »  0.25 
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/>  (0.75,0.25:0.75) 


0.55620 


P2(a1#«2:t)  -  P2 (0. 5, 0. 5: 0. 5)  =  0.39347 

Prob iem  (13 ) 

Gi  v  en 

3  1 

a  »  60°,  cr1  =  1,  cr2  =  1,  62  *  r  *  1 

“  b 

Find 

P  2  ( a  l »  a  2 ;  t  > 


This  problem  uses  the  same  elliptical  probability  distribution  and  an 
elliptical  area  over  which  the  elliptical  probability  is  to  be  integrated 
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thul  is  the  sume  size  and  shape  us  the  oilipticul  area  in  the  preceding 
prohlcm,  except  thut  lor  this  problem  the  elliptical  area  lias  been  rotated 
through  90°.  The  relationship  of  this  ellipse  to  the  probability  distribu¬ 
tion  can  be  seen  in  Fig,  C-l . 

C  ■  byv\x  +  b2crj2  -  (2 )  (3/2)  +  (2/3 )  (1/2 )  =  3(1/3) 

feig»i  (2X3/2) 

fl‘  C  =  3(1/3) 

b2crU  (2/3)  (1/2) 

°2  *  C  3(1/3) 

t  -  ±  »  — L_  =  0.3 

C  3(1/3) 

Pi'(al,a2\t )  =  (0.  9,  0. 1;  0. 3 )  «  0.34945 

Prob  lem  (14) 

Given  two  lines  of  position  which  intersect  at  an  angle  of  65°  and  which 
have  standard  deviations  of  700  1't  and  900  it,  respectively.  Find  the 
probability  that  the  indicated  position  is  within  a  quarter  of  a  nautical 
mile  of  the  true  position. 

by  the  given  quantities,  a  =  65°,  o-j  *»  7  00  ft.,  <?2  -  900  ft.,  bj  =  1, 

b2  =  1,  r  -  1520  ft. 


1  4(sin2  a)a?crj 

-  —  .  1  +  1 - s — t 

2  ,  n  n  .  n 


(cr2  +  O’2  )2 


1  4( 0. 82139 )  (TOO)2 (900) 2 

-  i  +  1 - 

2  (7002  +  9002 ) 2 


0.739 


0.261 


r2  ( 1 520) 2 ( 0. 82139 )  i 

cr\  +  o\  (  700 ) 2  +  (900 ) 2 


58  1,460 


PjUj,  a2;  t )  -  P2(0. 739,  (.261;  1.460)  =  0.77624 
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Prob I (  1 5 ) 

Given  Liu*  sunn'  conditions  in  the  preceding  problem,  but 
that  the  true  position  is  within  half  a  nautical  mile  i 

Qj  =  0.739 

u2  =  0,261 

r2  sin2  a  (3040 ) 2 (0.82139) 

t  =  — ~  —  -  ■■  "  =  —  . . — ~  =  5. 

cr\  +  cr\  (7  00)  2  +  (900) 2 

P 2 ( j . a 2 ; t )  =  P,(0. 739.0.261:5. 84)  =  0. 


the  probability 
s  desired. 

84 

99404 
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APPENDIX  D 


DERIVATION  OF  METHOD  1  FORMULAS 

1.  INTRODUCTION 

Method  1,  as  described  in  the  body  of  the  memorandum,  initially 
trouts  the  special  case  where  CTj  11  a2 .  Early  in  the  performance  of  this 
project  the  probability  curves  for  this  case  shown  in  Fig.  7,  Section  2, 
were  calculated  by  u  digital  computer  program.  However,  because  of  their 
limited  application,  little  use  could  be  made  of  them  until  the  further 
operations  involving  the  fictitious  cr*  and  a*  were  developed.  A  step  in 
the  development  of  these  functions  was  the  derivation  of  an  alternative 
set  of  formulus  for  a *  and  cr .  This  alternative  set  involves  the  computa¬ 
tion  of  auxiliary  functions  which  were  then  found  to  lead  to  the  formulas 
for  cr*  and  a*. 

This  appendix  discusses  in  order  the  development  of  the  probability 
curves  P(R/a,a)\  the  development  of  the  alternative  formulas  for  cr  and 
cry ;  and  then  gives  a  derivation  for  the  formulus  of  Figs.  9  and  10  of 
Section  2  that  utilizes  relationships  developed  in  Appendix  13  as  starting 
points. 

2.  DEVELOPMENT  OF  P(P/cr,a)  CURVES 

Given,  as  in  Appendix  B,  the  intersection  of  two  lines  of  position 
ut  an  ungle  ot.  however,  this  derivation  is  limited  to  the  special  case 
where  the  two  associated  standard  deviations,  cr1  and  a 2  ure  equal.  In 
this  special  case.  Equation  (5)  of  Appendix  B  may  be  simplified  to 

P  -  — L-  [[  e-U/2a*)(.*+.J)  dudv  .  (1) 

27 tct2  ■ 

Since  we  want  to  work  in  orthogonal  coordinates  we  sliull  proceed  to 
transform  to  new  orthogonal  coordinates  g  and  77  where  the  r\  axis  bisects 
the  angle  a  ( Fi g.  D- 1 ) .  The  choice  of  this  position  is  to  reduce  the  ulgcbra 
later  on. 
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ta-*  2 1  e  r-i  *  i 


FIG.D-1  TRANSFORMATION  OF  AXES 

Take  an  arbitrary  point  Q.  We  want  to  find  its  position  along  the 
x  and  y  axis  in  terms  of  £  and  T).  Take  the  vector  m  from  the  origin  0 
to  Q.  m  +  T)\^  where  l ^  and  l  are  unit  vectors  along  the  £  and  T} 

axes  respectively, 

/ t  and  l  are  unit  vectors  along  the  u  and  v  axes  respectively.  To 
find  the  projections  of  m  on  the  u  and  v  axes,  we  dot  it  with  the  respec 
tive  unit  vectors. 

-  -  -  -  —  ot  a 

u  8  m  '  lx  »  tjlg  '  lx  +  Tjl^  *  l  x  8  £  sin  —  v  tj  cos~ 

(2) 

“  '  Ty  +  Tflv  •  lf  =  -  sin  --+  v  cos~ 

(3) 

us  +  ti2  8  sin2  ~  +  17s  cosJ“"j  (4) 


dp 

2 //a- 

We  are 

now  in  I'l'RH 

1  a  r  0  r 1 

i  n  tii  0 

r  i  re  1  0  0  f  r 

ml  ins  1 

c.f 

V-*2 

I 


£  2  »  i  n2  -jr  +  r?  c  (i»2  j- 


i  n  adSjdri 


(6) 


HOMu!  UXOH  u  n  J  V  To  find  the  probabilities 

i  ||(.L'LrlilU‘  lfy>  U*CT  l!lO  circle. 


JS“  0 


77=0 


2;/cr4 


l  L 


2  H  i  «1  2  f  +7)2OUt2  -f 


sin  iid^d'r) 


(7) 


let 


^  n  r  cos  01  sin  “ 

j  •  y 


(8) 


uud 


/•  sin  ol  cos 


(9) 


(^‘/7'  e  J  (t  t)  drd *  =  TTTtf  (/r^ 


ao) 


uihI 


-  «/ 


/  cos2 

4> 

s  i  n2 

0 

— — 

/  si,,2 

a 

0 

sin2 

a 

~o 

1 

(11) 


We  therefore  have 


<£=«/'  2 


J0«Q 


'  =  «/ 


r“  0 


7  2  1 

coi  4  “  1  n  $ 

""Ta  ,  2  cf 

Bin"  *n*  f 


2770- 1 


0  2  rdrd<b  •  02) 
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Performing  the  first  integration 


fW/t  *  ^44)  * 

.  p  *  1  -  —  I  17  '  \"u  7  *3/  ti0  .  (15) 

7 T 

Jo 

This  integral  cannot  be  evaluated  analytically.  A  solution  was  ob¬ 
tained  by  a  digital  computer  program  and  the  results  have  been  presented 
in  Fig.  7,  Section  2,  and  numerically  in  Table  l,  Section  2. 

3,  ALTERNATIVE  FORMULAS  FOR  crf  AND  cry 

Additional  formulas  for  cr  and  utilizing  auxiliary  functions  were 
also  developed  in  the  course  of  the  analytical  investigations  of  this  study. 
The  additional  formulas  are  given  below,  followed  by  the  derivation  that  is 
developed  using  the  methods  and  notation  of  vector  algebra. 


cr  si  n  2/3 


cr  sin  2/3 


/3  =  urctan  cr^/cr^ 


(16) 


(17  ) 
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loin  wli  it'll 


ain  2/5 


2v  i r:  2 


..l 


The  cool'd i nuU  HyatomH 


°l  +  a2 

exhibited  ill  l; i B ■  n‘2- 


Si  net* 


and 


V,  =  «  * 

V2  «  R  1  £2 

A  =  *1^1  +  *2^2  1 

y  =  *,  COS  0  +  *2  sin  £ 

v  .  *J  cos( a  +  0)  *  *2  Sin(a  t  0) 


(18) 
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Now  if  U,0)  are  a  pair  of  vectors  where  A  -  (Al,Az)B  =  (B[tB2)  note  tli" 
identity: 


U,x,  +  Bj *2 ) 2  +  U,*i  +  B2x2)2  =  Ul2*?  +  2(A  •  B)*,*2  +lfi!2je2 

*  I x^A  +  x2B| 2 


Let 


A  =  [cos  ^/cjg,  cos  («  +  fiO/o-J 
8  =  [sin  0/cr  .  sin  (a  +  d)/ah] 


Choose  f?  so  that  4  '  B  B  0,  that  is. 


cos  fc"  sin  6  cos  (a  +  6)  sin  ( c*,  +  0) 

- — -  + - - - —————— 

2  t-A 


or  equivalently: 


•0-2  -  cr\ 

arctan  I -  tan  Otj  “(X 


to-2  +  cr? 


By  this  choice  of  0  with  the  help  of  (20): 

Ul 2  IbI 2  =  i A  x  fl|2  - 


Define 


=  cr  cos  /3 

i 

■h  *  cr  sin  (i 
Then  (21)  may  be  expressed 


—  {arctan  [cos  (2/5)  tan  ft]  "it} 

2 


Temporarily  define 


cr1  sin2  2/3 

H  =  cos  2 6  sin2  /3  +  cos  (2a  +  26)  cos2  /3 


(19) 


(20) 


(21) 


(22) 


(23) 


(24) 
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formulas  for  cr g  and  <Jf,  A  subsequent  derivation  that  starts  from  formulas 
of  Appendix  B,  however,  is  easier  to  follow. 

Formulas  (16)  and  (17)  of  this  section  may  be  shown  to  be  equivalent 
to  the  equations  developed  in  Appendix  B  (36)  and  (37).  To  do  so  it  is 
necessary  to  substitute  the  equivalents  for  cr  and  /3  given  in  this  ex¬ 
pression. 

The  first  step  is  to  square  Equation  (16)  to  eliminate  one  radial  in 
the  denominator  as  a  simplification 


sin2  2/3 


Now,  substituting 


2  1  -  /l  -  sin2  2/3  sin2  a 


cr  f  +  Crf 


sin2  2/3  = 


4crjcr2 


(o-2  +  cr2  ) ‘ 


2oial 


M  +  °i>J 


i  -i 


4  cr2cr2 


si  n2  OL 


2  \  2 


(cr2  +  a2) 


2al°\ 

(o-2  +  cr|  ) 


1  +  /I  - 


4o-2cr2 


(^2  +  a2)2 


sin2  a 


4crf  cr 


2  rr'l 


1  “ 


1  2 


(cr2  +!T2 


sin  Ot 


1  + 


i 


(cr2  +  a2)  X 


erf  +  Cri 


+  cr\  )  -  4  sin2  acr2cr2 


2 ( cr2  +  cr2  ) 2  sin2  a 
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} 


I  .1 


(<rj  +  er2  ) 


1  +  - - -  /(cr*  +  CTj  )  2  -  '4  Hin2  Owjcr2 

-  -  "?  +  ffl 


2  sin2  a 


(<r2  +  a2 )  t  /(crj  +  oj  )  2  ~  4  sin2  cr*cr 


2  sin  a 


QED 


Similarly,  except  for  changes  in  sign  we  may  show  that  cr2  from  this  section 
is  equivalent  to  Equation  (37)  of  Appendix  B. 


4.  DERIVATION  OF  FORMULAS  USED  FOR  CURVE  AND  NOMOGRAMS 
IN  SECTION  2 

In  the  body  of  the  memorandum,  Section  2,  these  formulas  were  given 
as  the  basis  for  Figs.  8  and  9.  This  curve  and  nomogram  permits  the  con¬ 
version  of  unequal  cTj  and  o'2  associated  with  the  intersection  angle  a  to 
a  fictitious  pair  of  equal  standard  deviations  cr*  and  a  fictitious  inter¬ 
section  angle  a*.  Following  these  conversions,  the  probability  curves  of 
Fig.  7,  Section  2,  may  be  entered. 


fl/cr 


a*  =  urcsin  (sin  2/3  sin  a) 

fl/2 


v/ff2  +  cr2 


c  sc  2/3 


The  derivation  starts  from  Equations  (34)  and  (35)  of  Appendix  B. 
In  any  error  ellipse  we  have 


and 


0-2cr2 
1  > 


<y 1  +  cri 


c-y2 


sin*  a 


°\  +  cr2 


sin  a 


Appendix  B  (34) 


Appendix  B  (35) 


(27) 


(28) 


:  i 
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We  assume  that  the  same  distribution  may  also  be  described  by 
fictitious  functions  cr* ,  a*,  and  OL* ,  Thus  substituting  in  the  formulas 
above 

cr*  4 

-  (29) 


ar*cr* 
*  y 


and 


r2  +  a2  - 


sin'  a* 


2  cr*2 
sin2  a.* 


(30) 


equating  the  right-hand  sides  of  Equations  (27)  with  (29)  and  (28)  with 
(30),,  we  obtain 


and 


o\a\ 


sin'  oi 


O’*  +  Oi 


r*  < 


sin'  a* 


2cr* 


sin2  a,  sin2  a* 

From  the  second  of  these  we  get 

sin2  fit  *  (cr?  +  cr? ) 

>  ,  l  i 

.  cr*  '  =  . . . . .  .  - - 

2  sin2  a 

Substituting  in  the  first  we  get 


(31) 


(32) 


(33) 


cr2cr2 


sin4  a*(<r2  +  cr2)2 


sin2  a  4  sin4  a  sin2  cl* 


sin2  a* 


sin  oi*  * 


4  sin4  occr2cr2 
sin2  a  ( cr2  +  o-2  ;  2 
20-^2  sin  01 
cr2  +  cr2 


(34) 


(35) 


(36) 


From  Section  3  of  this  appendix  we  obtain  the  following: 


Define 


arctan  crj/cr t 
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(37) 


i 

L 


( 


(38) 
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COMBINATION  OF  MULTIPLE  ELLIPSES 


Preceding  page  blank 
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APPENDIX  E 


COMBINATION  OF  MULTIPLE  ELLIPSES 


The  basic  elements  of  the  combination  of  several  ellipses  with  ran¬ 
dom  orientations  of  their  axes  have  been  given  in  Section  3  of  the  basic 
Memorandum  (see  Fig.  15).  Each  ellipse  ia  expressed  by  its  individual 
values  of  and  t r ^  and  the  orientation  angle  8  with  respect  to  the  arbi¬ 
trarily  selected  u>  and  z  axes.  The  first  step  in  the  combination  of  sev¬ 
eral  ellipses  is  to  transform  the  individual  values  of  cr *  and  a  to  values 
along  the  tv  and  z  axes.  Such  a  transformation  will  be  shown  to  involve 
also  a  cross-product  function  p  in  addition  to  the  new  standard  deviations 
cr  w  and  cr^  The  necessary  transformation  equations  are  the  standard  ones 
for  the  rotation  of  coordinate  axes. 

tv  *  x  cos  6  -  y  sin  9 

z  «  y  cos  6  +  x  sin  9 

or  equivalently, 

x  ■  u  cos  6  +  z  sin  6 

y  •  z  cos  6  -  w  sin  8  • 

The  original  distribution  in  terms  of  x  and  y  is  given  by 


at 


1 


2  vu  cr 


y 


Substituting  the  transformed  values  for  x  and  y,  we  obtain 


f(v, z ) 


-1  [lum  g  >i  LU  fl>2  + 

e  2  L 


(»  cos  8  -  »  iin  0) 
_2 
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Precedng  page  blank 


Considering  now  only  the  exponent  of  the  general  expression  of  an  ellip¬ 
tical  bivariate  distribution,  which  is 

_  1 
2(1  -  p 2) 


one  may  equate  these  two  exponents  and  solve  for  cr^,  crj  and  p.  First, 
multiply  out  the  trigonometric  terms  in  the  first  exponent 


1  fw2  cos2  9  2m  sin  9  cos  9  z2  sin2  9  z2  cos2  9  2wi  sin  9  cos  9  w2  sin2  o\ 
- -  + - - + -  +  - - + - 


rearranging  terms 


The  bracketed  coefficients  in  each  term  may  be  replaced  by  <4,  C,  and  B, 
respectively,  to  give  a  simplified  expression 


-  U24  +  2C»z  +  z  2fi] 
2 


This  simplified  expression  may  now  be  equated  to  the  general  expression 
of  the  exponent  and  solved  for  the  desired  functions  by  comparison  of 
similar  parts. 


i  i 


w2/l  -  2Cwz  +  z  iB 


1 


(1  “  P J) 


w  2  A  * 


cr 


2 

V 


1 

(1  -  p2)A 
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V 


\ 


similarly 


1 


(1  -  p2)B 


-2Cwi 


2  put 


o-  cr 

W  t 


Substituting  the  values  obtained  above  for  cr  w  and  cr 


p(  1  -  p 2)  ^AB 
(1  -  p2) 

C 

P  B  vtffl 

Thus  the  formulas  for  cr ,  cr  and  p  indicated  in  Section  3  have  been  de¬ 
rived.  The  subscript  t  used  in  Section  3  indicates  the  particular  ellipse 
by  number  where  t  ■  1,  2,  3  ...  n. 


Following  the  computation  of  the  functions  pertaining  to  each  sepa¬ 
rate  ellipse  with  respect  to  the  n  and  z  axes,  it  is  then  necessary  to 
combine  these  to  obtain  the  parameters  of  the  final  combined  error  ellipse. 
Since  the  final  error  is  the  sum  of  individual  errors,  then  the  moment- 
generating  function  of  the  final  distribution  is  the  product  of  the 
moment-generating  functions  of  the  individual  distributions  (since  the 
distributions  are  independent).  The  moment-generating  function  for  w  . , 


I  +  2p  |  (  cr  or  +  l2*2  \ 

•  -  e  2  '  1  *  1  *  *«/ 


Now 


nltj, t2) 


n  m((t 

i-1 


L  {l?(cr2  +  cr2  +  ...  cr2  ^  +  2t,t.(p.cr  cr  +  ...  p  cr  cr  ^ 

t  )  .  .  2  V\  "1  *2  *2 )  1  A1  *1  ‘1  ”  %  V 

l,r2'  e 
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where  m  (t  |,  1 2)  is  the  moment-generating  function  of  the  final  distribution 
and 


m(t  j , t  j )  =  e 


*  ipf‘ l**"./'.,  + 


Identifying  similar  parts,  we  have 


cr2  =»  1  (71 

1  i-\ 


2  a2 


1  ; 

~  1  pp*?* 


<y  cr  , 

i  w  i*  1 


Having  thus  obtained  the  parameters  of  the  final  combined  error 
ellipse  with  respect  to  the  arbitrary  w  and  2  axes,  it  is  further  desirable 
to  convert  to  the  x  and  y  axes  of  the  final  ellipse  to  remove  the  cross* 
product  function  p  ^  Formulas  for  this  transformation  may  be  found  in 
Ha  Id,  Ref.  8. 


cr  cr  ■  cr  cr 

•f  »/  */  */ 


a2  +  cr 2  »  cr  a 

X1  yf 


(Hald  19.8.10) 


(Ha  Id  19.8.11) 


These  equations  are  identical  in  form  with  (34)  and  (35)  of  Appendix  2. 
Thus  the  solutions  may  be  given  by  analogy  as 

cr2  «  —  [cr2  +CT2  +  /[cr2 +  cr2  \  -  4cr  CT ( 1  “  p2  )1 

*/  2L  ■/  */  /  \  */  ‘{I  */  '/  J 

cr2  »  -|cr2  +  cr2  -  /  (cr2  +  Cr2  \ “  4 <T  cr ( 1  "  O2 )  1 

»/  2L  •/  */  /  v  ■/  ‘fl  ■/  */  ■  J 

The  relationship  of  the  x  axis  of  the  final  ellipse  to  the  w  and  z  axes 
is  given  by 


tan  2 8  , 


2P/0\  cr 
/  w  /  / 


cr2  -  cr2 

*/  */ 


(Hald  19.8.7) 
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